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THE STABILITY OF ECCENTRICALLY 
STIFFENED CIRCULAR CYLINDERS 

VOLUME IV 

GENERAL INSTABILITY OF CYLINDERS 
HAVING LONGITUDINAL AND CIRCUMFERENTIAL 
STIFFENERS; AXIAL COMPRESSION 

By 

G. W. Smith, E. E. Spier, and E. A. Muscha 
General Dynamics Convair Division 
San Diego, California 

ABSTRACT 

This is the fourth of six volumes, each bearing the same report 
number, but dealing with separate problem areas concerning the stability 
of eccentrically stiffened circular cylinders* The complete set of documents 
was prepared under NASA Contract NAS8-11181. This particular volume deals 
with the general instability of simply supported, axially compressed circular 
cylinders which have longitudinal and circumferential stiffening (stringers 
and rings)* Analysis methods are presented in the forms of curves and a 
digital computer program. Since the contents of this volume are based 
upon a Donnell-type small-deflection theory, the proposed methods should be 
used in conjunction with empirical knock-down factors to account for the 
effects of initial imperfections. In addition, the Donnell assumptions pre- 
clude application to non-axi symmetric buckle patterns where the number of 
circumferential waves is small. 
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r+s 


A 

s 


A X1’ A 22’ 

A A 
12’ 33 



R 

F 



DEFINITION OF SYMBOLS 
Definition 

Cross-sectional area of single ring (no 
cylindrical skin included) ; (see Table III 
and its notes) . 

Cross-sectional area of single ring augmented 
by an effective width of skin which is considered 
to essentially behave as part of the ring (see 
Table III and its notes)* 

Cross-sectional area of single stringer (no 
cylindrical skin included) ; (see Table III and 
its notes)* 

Elastic constants (see Table III). 

Ring spacing. 

See equations (2-14) and inequality (2-15) . 
Stringer spacing. 

Thickness of integral longitudinal stiffener 
(see Table III). 

Correction factor defined by equation (3-5)* 
Experimentally determined constant [see equation 
( 1 - 1 )]. 

Quantity defined by equation (2-25). 
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DEFINITION OF SYMBOLS 
(Continued) 


Symbol 

C hy 

V 


py 

C 11’ C 22’ C 12 ,C 21 


c 

D 


D. 


D 11’ D 22’ D 12’ D 33 

E 

E F 

G 

G 


Definition 

Quantity defined by equation (2-29). 

Quantity defined by equation (2-26). 

Quantity defined by equation (2-30). 

Eccentricity coupling constants (see Table III). 
Thickness of sandwich core. 

Diameter of middle surface of basic cylindrical 
skin in conventional skin-stringer-ring con- 
structions; Parameter defined in equations (A-l) . 
Transverse shear stiffness of sandwich wall 
[see equations (A-l)]. 

Elastic constants (see Table III). 

Young's modulus. 

Young's modulus of sandwich faces. 

Modulus of elasticity in shear. 

Shear modulus of sandwich core. 

Distance between middle surfaces of faces in 
sandwich cylinder. 

Depth of integral longitudinal stiffener 
(see Table III). 

See Figure 2. 

See Figure 3. 
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Symbol 

I 

I 

r 

I 

r+s 


I 

x 


I 

y 


i 



L 

M 

m 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 
Moment of inertia. 

Centroidal moment of inertia for single ring 
(no cylindrical skin included); (see notes 
following Table III). 

Centroidal moment of inertia of single ring 
augmented by an effective width of skin which 
is considered to essentially behave as part of the 
ring (see Table III and its notes). 

Running centroidal moment of inertia of effective 
shell wall cross section lying in plane normal 
to axis of revolution (see Table III and its notes)'. 
Running centroidal moment of inertia of effective 

shell wall cross section lying in radial plane. 

• "til 

= V -1 ; symbol denoting the i quantity of an 

array of values. 

Parameter defined by equation (3-7) . 

Sandwich loading parameter defined in equations 
(A-l). 

Overall length of cylinder (see Figure l). 

Overall bending moment. 

Number of axial half-waves in buckle pattern. 
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Symbol 

N 

s 


n thiel 


(**THIEL 

N 

x 


c 



n 


R 



DEFINITION OF SYMBOLS 
(Continued) 

Definition 


Number of discrete load points for the case 
shown in Figure 4(a); (Number of stringers). 

Loading parameter defined in equations (2-3)) 
(positive for tensile loading). 

= - N thIEL’ (P os **i ve * or compressive loading) . 
Applied longitudinal tensile running load acting 
at the centroid of the effective skin-stringer 
combination. 

Applied longitudinal compressive running load 
acting at the centroid of the effective skin- 
stringer combination ( = - N ) . 

Number of circumferential full-waves in buckle 
pattern. 

Radius to middle surface of basic cylindrical 

skin in conventional skin-stringer-ring constructions 

Radius to centroid of ring depicted in Figure 4; 

Mean radius of sandwich cylinder; Mean radius of 
isotropic cylinder. 

Parameter^ defined in equations (A-l). 

Thickness of basic cylindrical skin. 

Thickness of single face of sandwich wall. 

Skin thickness of corrugated wall. 
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Symbol 


t 

x 



U 

u,v,w, 

W 


x,y,z 

z 



P 



DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Thickness of appropriate sraeared-out area of 
cross section lying in plane normal to axis 
of revolution (see Table III and its notes). 
Thickness of appropriate smeared-out area of 
cross section lying in radial plane (see Table 
III and its notes). 

Parameter defined by equation (2-10). 
Reference-surface displacements (see Figure l). 
Force depicted in Figure 4(a). 

Running load depicted in Figure 4(b) • 
Coordinates (see Figure l). 

Parameter defined by equation (2-2). 

Parameter defined in equations (A-l). 

Parameter defined by equation (2-7). 
Eccentricity (see Figure 2 and Table III). 
Eccentricity (see Figure 3 and Table III). 
Parameter defined in equations (2-3). 

Parameter defined in equations (2-3). 

Parameter defined in equations (2-3). 

Radial deflection for the points of load 
application shown in Figure 4(a). 
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DEFINITION OF SYMBOLS 
(Continued) 

Definition 

Radial deflection due to uniformly distributed 
running load shown in Figure 4(b). 

Parameter defined in equations (2-3). 

Parameter defined in equations (2-3)<t 
Half-angle between discrete load points depicted 
in Figure 4(a). 

Poisson's ratio. 

Poisson's ratio for faces of sandwich cylinder. 
Local centroidal radius of gyration for effective 
shell wall cross section lying in a plane which 
is normal to the axis of revolution (see Table 
III and its notes). 

Local centroidal radius of gyration for effective 
shell wall cross section lying in a radial plane 
(see Table III and its notes). 

Parameter defined in equations (2-13). 

Total peripheral length of corrugation center- 
line. 

Critical compressive buckling stress of isotropic 
cylinder. 

Compressive stress obtained by dividing the 
longitudinal compressive running load by t^ 
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SECTION 1 
INTRODUCTION 

The contents of this volume deal with the general instability (see 
GLOSSARY, Volume I Cl]) of simply supported, axially compressed circular 
cylinders having both stringers and rings. Designing for the prevention 
of this mode of instability usually centers around the choice of a suitable 
criterion to establish dimensions for the circumferential stiffeners. In 
the past, a number of rather crude empirical formulas have been proposed 
for this purpose. One of the earliest of these was proposed by the 
Guggenheim Aeronautical Laboratory of the California Institute of Technology 
(GALCIT) as an outgrowth of their tests on small-scale cylinders C 2] . 

Shortly thereafter, the Polytechnic Institute of Brooklyn Aeronautical 
Laboratory (PIBAL) proposed a different criterion based on their own test 
results from similar specimens [ 3 ]. Shanley L4-] then drew upon both the 
GALCIT and the PIBAL data to generate the following empirical formula for 
the minimum ring stiffness required to prevent general instability in 
stiffened cylinders subjected to pure bending: 

C.MD 2 

El = — (1-1) 

r a 

where 

E = Young's modulus 

1^ = Centroidal moment of inertia for ring 

C^, = Experimentally determined constant 

M = Overall bending moment 

D = Cylinder diameter 

a = Ring Spacing 


Numbers in brackets t ] in the text denote references listed in SECTION 7. 
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Although derived specifically for the case of pure bending} this formula 
has been widely used for cases of axial compression by considering the 
peak running load intensity to be the controlling factor. Still another 
criterion was suggested by Becker C 5l in 1958. This approach employs 
certain geometrical features of the inward-bulge along with an estimate 
of the elastic restraints afforded by the rings. In general} all of these 
approaches constitute oversimplifications of the problem in that they do 
not recognize all the important variables involved. Engineers have long 
been wary of these criteria and have hedged their ring designs through the 
use of generous safety factors and extensive proof-testing. 

In recent years it has become increasingly popular to attack this 
problem by means of orthotropic shell theory. For example, applications 
have been made [ 6*1 of the small-deflection formulation presented by 
Thielemann in reference 7. The methods presented in the sections to follow 
are based on an extension of this approach to include the effects of finite 
cylinder length as well as stringer and ring eccentricities. These in- 
fluences have also been included in an alternative formulation developed by 
Block, Card, and Mikulas in reference 8. Furthermore, Block [9] has recently 
published a paper which treats the buckling of eccentrically stiffened 
perfect orthotropic cylinders in which prebuckling deformations, load 
eccentricity, and ring discreteness are all considered. That paper points 
the way for further improvement of the methods presented in the sections to 
follow. It should be noted here that, under NASA Contract NAS 8-9500, the 
Lockheed Missiles and Space Company has likewise published a solution ClO] 
which accounts for prebuckling deformations and load eccentricities (end 
moments). However, their development is limited to configurations which 
are stiffened only in the longitudinal direction (no intermediate rings). 

The analysis methods presented in this volume are given in the forms 
of curves and a digital computer program. In order to keep the number of 
curves within reasonable bounds, they embody a number of simplifications 
which result is some loss of rigor. Therefore, these curves are to be used 
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primarily for the purposes of preliminary sizing, rough checking, and the 
study of trends. To meet the more stringent requirements of a final 
design, one should employ the digital computer program of SECTION 6. 

This program can be used to obtain single-point solutions and/or additional 
curves which more closely apply to the selected configurations of interest. 

It should be kept in mind that this volume deals only with the general 
instability mode of failure. Other possible modes must be separately checked*, 
The most important of these other possibilities is the panel instability 
(see GLOSSARY, Volume I Cl]) mode which is treated in Volume III [ll] . 

The theoretical foundation for the methods given here lies in small- 
deflection orthotropic shell theory. Therefore, it is recommended that 
these methods be used in conjunction with empirical knock-down factors in 
order to account for the influences of initial imperfections. Volume V 
[l2l presents a practical criterion in this regard. Since most practical 
stiffened cylinders are ’’effectively thick", the related reductions will 
generally not be nearly so severe as those encountered for thin-walled 
isotropic cylinders. 

It is also important to note that the basic buckling equation of this 
volume is based upon Donnell-type simplifications [l3]. As a result, the 
methods given here cannot be applied when the instability manifests itself 
in a non-axi symmetric buckle pattern having a small number of circumferential 
waves. The rule-of-thumb guideline is offered here that these methods 
should be considered inapplicable for cases where 

0 < n < 2 (1-2) 

Throughout this volume, only the case of classical simple support 
is considered. In addition, it is always assumed that the behavior is 
elastic. To consider critical stresses which lie above the proportional 
limit, one must employ iterative procedures which recognize the shape of 
the applicable stress-strain curve. 
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SECTION 2 
EQUATIONS 


2.1 GENERAL 

The following orthotropic cylinder equation provides the basis for 
the methods given in this volume: 


(^THIEL^ 


1 + V7 P 2 + r P 4 a P 4 (z ) 2 


4oP 4 [l + 2T) s 0 2 + P 4 ] 


( 2 - 1 ) 


where , 


Z = 


1 - 


C + C 
11 22 


12 


2a(A D ) 4 ^ 2 6 2 
2222 ; p 


2 “ A 22 (D 22 /A U ) 


1/2 


21 


2<X ^ A 11 I) 22^ ^ 


( 2 - 2 ) 


A detailed derivation of these relationships is given in reference 14 where 
the coordinate system shown in Figure 1 was used. Some general background 
information concerning equations (2-l) and (2-2) is given in Volume I [ l] . 
As noted there, these equations have been written in rather compact, in- 
structive forms through the introduction of the following parameters, most 
of which were first proposed by Thielemann [ 7 ]: 


/s 

l THIEL 



N R 

N - _ -2L_ 

THIEL _ 2 




V A H A 22 


(2-3) 
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( P 12 + 2 P 33 ) 
^ D 11 D 22 


Y 


P 11 A 11 

° 22 A 22 



a 


2Rm 2 u 2 A 




1/2 


(2-3) 
cont 1 d 
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The various A. 's, D. . 's, and C. . *s of equations (2-2) and (2-3) are 
very important fundamental constants. The physical significance of these 
constants is discussed in Volume I Cl]. The A. ,'s and D. . r s are usually 

1 j 

referred to as elastic constants while the C. .'s might be identified as 

^ V 

eccentricity coupling constants. 


2.2 AXI SYMMETRIC MODE 

In order to apply equations (2-1 ) and (2-2) to practical structures, 
both axisymmetric and non-axisymmetric (checkerboard) buckling modes must 
be considered. Specialization to the axisymmetric case is achieved by 
allowing {3 to increase without bound (p -* "). This yields the following 
result: „ 

r 

C 

1 - 


im {N, 
B -» » V 


Lim 

P 


THIEL/ = 4a 


_ JL + a 


/c 


12 


2aA 22 (D 22^ A ll* 


1/2 


(2-4) 


By substituting equations (2-3) into (2-4), one can then obtain the following 
equation for axisymmetric buckling: 


M . • ' 


2 2 n 
it D 


11 


2 2 d 2 * 
n t R A 


22 


2 2 

1 - C 12 R 


(2-5) 


It is easily shown that, for m > 0 , the second derivative of (N x ) 


respect to 


m 


is always positive. Hence, a relative minimum for 


with 


(N ) can be located by differentiating equation (2-5) with respect to m 
x c 

and equating the result to zero. It is found that the derivative vanishes 


when 


where 


m = 


■ H 


1/4 


( 2 - 6 ) 


V 


r 2 a 


22 


^ D ll + 


C 12 2 / A 22) 


(2-7) 
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The quantity will always be a positive number. The four roots to 
equation (2-6) can therefore be expressed as follows: 



For applications to actual structures, only the single, real, positive 
root, m^, is of interest. For cases where m^ i 1, the critical loading 
for axi symmetric buckling can be found by making the substitution m = 
in equation (2-5). The condition m < 1 is a physical impossibility for 
the problem under discussion. Therefore, whenever < 1, the critical 

loading for axisymmetric buckling is found by making the substitution m = 1 
in equation (2-5). 

Detailed derivations of the foregoing relationships are given in 
reference 15 which also presents the background information concerning the 
digital computer program of SECTION 6. This program makes use of the above 
equations in the computation of output listings and the generation of 
automatically plotted critical strains. 

2.3 NON- AXISYMMETRIC MODE 

To obtain the critical loading for non-axi symmetric buckling, equation 
(2-l) must be minimized with respect to two wave-type parameters. In this 
volume, a and p were selected for this purpose although m and n would 
have been equally suitable. Therefore, the partial derivative of equation 
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(2-1 ) was taken with respect to a and the result was set equal to zero. 
From this it was found that the derivative will vanish whenever 


a = 


< 


l + 2T] p Vrp 2 + tE 4 ][i + 2T1 s p 2 + p 4 




(2-9) 


where 


0 _ ( c n ^ c ^) 

u — 1 yy 


ril * C 22 / ^12 "21 

2(A D 2A (D /A 2(A D 

22 22 ; P ^22 V 22'li ; 1 11 22 ; 


( 2 - 10 ) 


22 ' 22 ' 11 ' 


Here again, only the single, positive root is of interest to the problem 
under discussion. No theoretical consideration was given to the question 
of insuring that this a corresponds to a relative minimum rather than a 
relative maximum or an inflection point. Such safeguards were incorporated 
in the numerical procedures into which these relationships were injected. 

In addition, the minimization with respect to p is accomplished solely by 
numerical procedures. The digital computer program of SECTION 6 accomplishes 
this by screening over practical p values in small increments of P , always 
retaining the lowermost loading encountered. Constraints are imposed which 
disallow buckle patterns for which m < 1 and/or n is greater than zero 
but less than a set value selected by the analyst (usually unity ) * 

A detailed derivation of equation (2-9) is given in reference 15 which, 
as already noted, also presents background information concerning the digital 
computer program of SECTION 6. 

2.4 PREFERRED MATHEMATICAL FORMS 

In order to better understand the basic equations of this volume, and 
to facilitate the presentation of design data in the most practical manner, 
it is helpful to rewrite the equations in different forms from those 
presented in the foregoing. Reference 15 gives all the detailed mathematical 
steps involved in these transformations. Only the major results are given 
here. 
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2.4.1 GENERAL 

One type of useful reformulation of the basic buckling equation 
(2-l) can be achieved through the substitution of equations (2-3) and 
subsequent algebraic manipulation. This leads to equation (2-11 ) which 
provides the basis for the STIFF option computations of digital computer 
program 4267 (see SECTION 6). 

To obtain the additional reformulation expressed as equation (2-12) , 
it is helpful to define the following quantities: 




11 A 11 


y 




22 A 22 


(2-13) 




I 


*33 lD 12 * “33 ' 


Recognizing the physical basis [ 1 ] for the elastic constants D. . and A. , 

x* iJ ^ 3 

it is clear that the first two of equations (2-13) are quite logically 
selected as effective local radii of gyration in the axial and circumferential 
directions. The last of these three equations is an artif icial^contrived 
definition which has been made in the interest of convenience* It should 
be noted, however, that this equation does display some similaritv-of-f orm 
to the equations for p and p . It was this likeness which led to the 

y 

use here of the notation p * Additional clarification of the terms in- 

xy 

volved in equations (2-13) can be obtained from a study of Table III and its 
related notes. 

In addition to the foregoing, equation (2-12) is also based on the use 
of the following simplified formulas for the elastic constants of interest: 
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Table III and its notes furnish helpful clarifying information concerning the 
quantities which appear in these expressions. Permissible values for B A12 
are bounded as follows: 

0 s B A12 * 1 (2-15) 

As noted in reference 15, little error is introduced if this quantity is 
taken equal to 0 o 50 for all cases. Therefore, this value was selected for 
use in the present volume. This practice, together with the substitution 
of equations (2-13) and (2-14) into equation (2-11), permits one to rewrite 
the basic buckling equation in the instructive form of equation (2-12). 

This expression provides the basis for the RATIO option of digital computer 
program 4267 (see SECTION 6). 

2.4.2 AXI SYMMETRIC MODE 

The procedures described in SECTION 2.4.1 were likewise applied to 
rewrite equations (2-5), (2-6), and (2-7) in alternative, preferred 
mathematical forms. The following expressions were obtained for use in 
the RATIO option computations of digital computer program 4267 (see 
SECTION 6): 
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Since equations (2-5), (2-6 ) ? and (2-7) are already in a suitable form 
for the STIFF option computations, no additional transformations of 
these equations are required. 

2.4.3 NON-AXISYMMETRIC MODE 

The procedures described in SECTION 2.4.1 were likewise applied to 
rewrite equations (2-9) and (2-10) in alternative, preferred mathematical 
forms. Equations (2-18) and (2-19) were obtained for use in the RATIO 
option computations of digital computer program 4267 (see SECTION 6). 
Equations (2-20) and (2-2l) were obtained for use in the STIFF option com- 
putations of the program. 

2.5 SPECIAL GEOMETRICAL CONSIDERATIONS 

The large number of geometrical values appearing in the foregoing 
equations makes it necessary to introduce a number of simplifying assumptions 
to arrive at a suitable format for general instability design curves. 

Several such simplifications have been selected based on practical con- 
siderations relative to the cross sections which lie in planes normal to 
the axis of revolution. A portion of such a section is depicted in 
Figure 2. It should be noted that b denotes the circumferential stringer 


2-9 

GENERAL DYNAMICS CONVAIR DIVISION 



GENERAL DYNAMICS CONVAIR DIVISION 



9oo-Z9-oaa-ociD 



L DYNAMICS CONVAIR DIVISION 



GDC- DDG- 6 7-006 



GDC- DDG- 6 7-006 




Figure 2 - Sample Wall Cross Section 
Lying in Plane Normal to 
Axis of Revolution 


spacing. In addition, the area of a single stringer (no skin included) 
will be identified by the symbol A . The eccentricity is represented as 
. By assuming that 

(a) All of the skin and stringer material is fully effective 
and 

(b) The stringer centroid (no skin included) lies at its mid- 
height, 

one can employ the simple laws of mechanics to arrive at the following 
formula for the absolute magnitude of the eccentricity: 



( 2 - 22 ) 
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where 


t = 
x 


A 

f ♦* 


Introduction of the assumption that h^ » (t/2) permits one to rewrite 
equation ( 2 - 22 ) as 


(I M - uM 


(2-23) 


(2-24) 


It then becomes convenient to define two new quantities, and Cp^, 

as follows: 


hx 


(2-25) 


px 


(2-26) 


The quantity is the local centroidal radius of gyration of the skin- 

stringer combination. For practical skin-stringer configurations, the value 

C = .33 will usually approximate the actual conditions reasonably well. 

P ^ 

Therefore, the buckling curves of SECTION 4 are based upon this value. 

These curves are also based upon the relationship 


hx 


= (i\ l(y*h il 

M 


(2-27) 
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which is obtained through the substitution of equation (2-25) into (2-24). 
The computations performed within digital computer program 4267 make use 
of the C^ x and values to arrive at reasonable measures of eccentricity- 
based on specified P values. That is, z is found from the following 
equation which results from the substitution of equation (2-26) into 
(2-25) : 



Although the buckling curves of SECTION 4 incorporate the several approx- 
imations indicated above, it should be noted that digital computer program 
4267 features several options which allow one to obtain improved accuracy 
for given configurations. For example, one may choose to input actual C^ x 

and C values for selected structures. Furthermore, an option is provided 
px 

which permits the insertion of primary input geometry in the form of elastic 
constants and eccentricity coupling constants. 

Consideration will now be given to cross sections which lie in radial 
planes. A portion of such a section is depicted in Figure 3. The symbol 
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a denotes the spacing between rings while the eccentricity value of 
interest is identified as z . Unlike the case for the section shown 

y 

in Figure 2, the usual proportions encountered in skin-stringer-ring 
constructions make it improper to assume that the complete skin-ring com- 
bimation is fully effective. Normally, only a relatively short length of 
skin can be considered to essentially behave as part of the ring cross section. 
Therefore, no relationship corresponding to equation (2-24) is involved in 
the ring-related simplifications embodied in the buckling curves of SECTION 4. 
These simplifications are introduced only through the terms C. and C 

r J 

which are defined as follows: 



(2-29) 



(2-30) 


The quantity p is the local centroidal radius of gyration of the ring- 

y 

effective skin combination. Substitution of equation (2-30) into (2-29) 
permits one to express the absolute magnitude of the hoop eccentricity as 
follows : 



(2-31) 


To generate the buckling curves of SECTION 4, this equation was used to 
arrive at reasonable measures of j z x j based upon the specified P^.. The values 


C hy = * 425 
C Q = .4 

py 


(2-32) 


were used for this purpose. These were selected as reasonable estimates for 
the proportions likely to be found in practical skin-stringer-ring con- 
structions. Note however that digital computer program 4267 features an 
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input format which allows one to obtain improved accuracy for given 
configurations. For example, one may choose to input actual and 
Cpy values for selected structures. Furthermore, as already noted, an 
option is provided which permits the insertion of primary input geometry 
in the form of elastic constants and eccentricity coupling constants. 
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SECTION 3 
ANALYSIS METHODS 


The methods given in this volume deal only with general instability 
(see GLOSSARY, Volume I Cl]), No consideration is given here to the panel 
instability mode (see GLOSSARY, Volume I Cl]). The methods employ the 
smearing-out technique whereby discrete stiffness values are averaged over 
the entire surface of the cylinder. One must therefore exercise engineering 
judgement to prevent misapplication to configurations having excessively 
large stiffener spacings. In addition, it should be noted that only classical 
theoretical solutions are employed. Influences from initial imperfections 
are completely ignored. It is therefore recommended that the values obtained 
from this volume be reduced in accordance with the knock-down criterion of 
Volume V [l2]. 

The methods presented here for the analysis of general instability 
basically consist of the following: 

(a) A collection of buckling curves (see SECTION 4) which embody a 
number of simplifications which are reasonable for conventional 
skin-stringer-ring configurations. 

(b) Digital computer program 4267 (see SECTION 6) which can be used 

to obtain single-point solutions and/or additional buckling curves. 
The input options of this program enable one to obtain particular 
solutions of greater accuracy than can be obtained from the curves 
of SECTION 4. These curves sacrifice accuracy in the interest of 
generality and practicality. The computer program also enables one 
to analyze corrugated and waffle configurations, neither of which 
are covered by the given buckling curves since all of these plots 
apply only to conventional skin-stringer-ring arrangements and 
proportions. 

The simplifying assumptions embodied in the buckling curves of SECTION 4 
are as follows: 

(a) The entire skin-stringer combination is assumed to be fully 
effective. 

(b) The stringer centroid (no skin included) is assumed to lie at 


its mid-height. 

(c> h x »! 
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(d) 

(e) 

(f) 

(ff) 

(h) 

(i) 


C Px 

sr 

.33 

C. 

— 

.425 

hy 

c 

rr 

.4 

py 

D _ 

— 

0 

33 

CM 

i— 1 

ft 

= 

°21 = ° 

C 12 

= 

C 21 = ° 


Assumptions (h) and (i) are justified on the basis that the constants D^, 
^21’ ^12’ ant * ^21 involve available mechanisms for interplay which are 
dependent upon Poisson' s-ratio influences. For usual stringer and ring 
arrangements, this interplay is of a minimal nature since it occurs primarily 
in the localized regions where the stiffeners intersect. However, when the 
stringer and ring spacings are very small (as for waffle configurations) 
and/or the contributions from stiffener rigidities approach that of the 
skin (also as for waffle configurations) , non-zero values should be used 
for D.^, ®21* ^12 an< * ^21* suc ^ cases, one must resort to point-solutions 

using digital computer program 4267, 

It is important to recall here that a Donnell-type theory furnishes the 
basis for the methods given in this volume. Therefore, these methods cannot 
be applied to cases of non-axisymmetric buckling where the number of 
circumferential full-waves is small. As a rule-of-thumb, it is suggested 
that the methods be considered inapplicable where 


0 < n < 2 (3-1) 

All of the curves given in SECTION 4 completely ignore this restriction. 
Therefore, to insure that this condition is not violated, one should supple- 
ment the plotted data with appropriate checks from digital computer runs 
using program 4267. However, most practical configurations likely to be 
encountered will not display buckle patterns with 0 < n < 2. Hence it 
should not prove necessary to make the suggested check for every configura- 
tion investigated. When a large number of candidate designs are to be 
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studied, it will usually be reasonable for one to assume that the wave- 
number restrictions are satisfied so that the foregoing check need only 
be made as a final operation for a few selected cases. 

It is further suggested that, in applying the methods of this volume, 
one give some consideration to the influences which shallow rings and/or 
wide stringer spacings might have on the elastic constant A 2 2* For the 
axisymmetric buckling mode it seems clear that the A value should reflect 
the difference in flexibility between the discretely loaded ring of Figure 4(a) 
and its uniformly loaded counterpart shown in Figure 4(b). For this purpose, 
a correction factor can be established which is simply a ratio of these 

two flexibilities. 
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The following formula C 16 3 can be used for the radial deflection of the 
discretely loaded ring: 


R 



r 


* 

WR 3 

1 / 

0 sin0 cos0 \ 

1 

2EI 

r 

sin 2 0 ' 

2 + 2 / 

~ 0 

- 


(3-2) 


where 

= Radial deflection at points of load application. 

0 = Half-angle between discrete load points, radians. 

The radial deflection for the uniformly loaded case is denoted 6^ and 
the following formula for this quantity is easily derived: 


w R 
c 

J R " A E 
r 


(3-3) 


To establish a basis for comparison between these two situations, the 

following relationship between W and w is employed: 

c 


w 


2ttR 
— r— w 
N c 
s 


(3-4) 


where 


N g = Number of discrete load points for the case of Figure 


4(a); (Number of stringers) 
The factor may then be defined as follows: 

6 , 


R 


Cr 4 R + 6 R 


a r 

T~ + 1 
6 R 


(3-5) 
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By substituting equations (3-2), (3-3), and (3-4) into equation (3-5), 
the following result is obtained: 


where 



(3-6) 


K 


nM 2 


sin 6 


sinQ cos8 
2 



1 

e 


(3-7) 


and 

9 = JT (3-8) 

s 

Figure 5 presents plots of C^ versus (R A^/I^) for selected constant 
values of N g . The form of equation (3-7) leads to the usual numerical 
difficulties that one encounters when handling small differences between 
large numbers. Therefore, the curves of Figure 5 were plotted from digital 
computer results obtained to double-precision accuracy. 

It is thought that the most logical use of the C^ value is as a 
correction factor to the quantity (l/A ^) appearing in any terms which do 
not vanish as P "* * . However, such corrections have not yet been in- 
corporated into either the buckling curves of SECTION 4 or digital computer 
program 4267. Therefore, Figure 5 can only be used here as a means to check 
the applicability of the given methods. It is suggested that, when C^ * 0.95, 
one might conclude that the discreteness mechanism depicted in Figure 4 is 
of negligible importance. On the other hand, lower values would indicate 
a need to apply engineering judgement in evaluating results obtained from 
the curves of SECTION 4 and digital computer program 4267. However, it will 
be found that, for most realistic skin- stringer-ring constructions, the 
geometric proportions will be such that C^ is essentially equal to unity. 
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SECTION 4 
BUCKLING CURVES 


The curves of this section present critical strain values for axially 
compressed circular cylinders having eccentric stringers and rings. To 
make proper use of these curves, one should refer to the instructions 
furnished in SECTION 3, "ANALYSIS METHODS". All of these curves were 
developed by using digital computer program 4267 (see SECTION 6) in 
conjunction with an automatic plotting machine. The machine located in- 
dividual points through which the curves were drawn by hand. Since the 
plotting machine does not have the capability to print out lower case 
letters, the symbols t, x, and the y of t are shown in upper case notation. 

y 

The following information is furnished to clarify the meanings of all 
terms appearing on the plots: 


T BAR X = t 

x 


T BAR Y = t 

y 


T = t 


Thickness of appropriate smeared-out area of 
cross sections lying in planes normal to the 
axis of revolution (see Table III and its notes). 
Thickness of appropriate smeared-out area of 
cross sections lying in radial planes (see 
Table III and its notes). 

Thickness of basic cylindrical skin in con- 
ventional skin-stringer-ring constructions. 

To facilitate application to other configura- 
tions and to cases where buckling of the iso- 
tropic skin panels precedes general instability, 
it is helpful to note that this quantity enters 
into the given curves through both the elastic 

constant A and equation (2-27) . 
oo 
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L s Overall length of the cylinder. (Note 

that this quantity is NOT the ring spacing). 

R = Radius to middle surface of basic cylindrical 
skin. 

p = Local circumferential radius of gyration. 

y 

This quantity may be computed from 



as applicable (see Table III and its notes). 


RHO X = = Local longitudinal radius of gyration. This 

quantity may be computed from the equation 

(4-2) 

(See Table III and its notes). 

SIGMA BAR = a = Critical compressive buckling stress obtained 

by dividing the critical longitudinal com- 
pressive running load by t x . 

E = Young's modulus. 

Table I lists the families of curves provided here for the case of external 
stringers and internal rings. Table II lists the families of curves 
provided here for the case of internal stringers and internal rings. 
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TABLE I - Table of 

Contents for Curves of 



Critical 

Axial Compression for 



General Instability 

in Cylinders 



With Eccentric Orthotropic Stiffening 



(External 

Stringers 

and Internal Rings) 


Figure 

/ 1 \ 
( x j 

( L \ 

K) 


Number 


U/ 


Page 

6(a) 

1.2 

1.0 

0.1 

4-5 

6(b) 

1.2 

1.0 

0.5 

4-6 

6(c) 

1.2 

1.0 

1.0 

4-7 

6(d) 

2.0 

1.0 

0.1 

4-8 

6(e) 

2.0 

1.0 

0.5 

4-9 

6(f) 

2.0 

1.0 

1.0 

4-10 

6(g) 

3.0 

1.0 

0.1 

4-11 

6(h) 

3.0 

1.0 

0.5 

4-12 

6 ( i ) 

3.0 

1.0 

1.0 

4-13 

6<j) 

1.2 

2.0 

0.1 

4-14 

6(k) 

1.2 

2.0 

0.5 

4-15 

6(1) 

1.2 

2.0 

1.0 

4-16 

6(m) 

2.0 

2.0 

0.1 

4-17 

6(n) 

2.0 

2.0 

0.5 

4-18 

6(o) 

2.0 

2.0 

1.0 

4-19 

6 (p) 

3.0 

2.0 

0.1 

4-20 

6 (q ) 

3.0 

2.0 

0.5 

4-21 

6 (r) 

3.0 

2.0 

1.0 

4-22 

6 (s) 

1.2 

3.0 

0.1 

4-23 

6 (t ) 

1.2 

3.0 

0.5 

4-24 

6 (u) 

1.2 

3.0 

1.0 

4-25 

6 (v) 

2.0 

3.0 

0.1 

4-26 

6(w) 

2.0 

3.0 

0.5 

4-27 

6 (x) 

2.0 

3.0 

1.0 

4-28 

6(y) 

3.0 

3.0 

0.1 

4-29 

6(z) 

3.0 

3.0 

0.5 

4-30 

6 (aa) 

3.0 

3.0 

1.0 

4-31 
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TABLE II - Table of Contents for Curves of Critical 
Axial Compression for General Instability 
in Cylinders With Eccentric Orthotropic 
Stiffening (Internal Stringers and 
Internal Rings) 


Figure 

Number 



7(a) 

1.2 

1.0 

0.1 

4-32 

7(b) 

1.2 

1.0 

0.5 

4-33 

7(c) 

1.2 

1.0 

1.0 

4-34 

7(d) 

2.0 

1.0 

0.1 

4-35 

7(e) 

2.0 

1.0 

0.5 

4-36 

7(f) 

2.0 

1.0 

1.0 

4-37 

7(g) 

3.0 

1.0 

0.1 

4-38 

7(h) 

3.0 

1.0 

0.5 

4-39 

7(i) 

3.0 

1.0 

1.0 

4-40 

7(j) 

1.2 

2.0 

0.1 

4-41 

7(k) 

1.2 

2.0 

0.5 

4-42 

7(1) 

1.2 

2.0 

1.0 

4-43 

7(m) 

2.0 

2.0 

0.1 

4-44 

7(n) 

2.0 

2.0 

0.5 

4-45 

7(o) 

2.0 

2.0 

1.0 

4-46 

7(p) 

3.0 

2.0 

0.1 

4-47 

7(q) 

3.0 

2.0 

0.5 

4-48 

7(r) 

3.0 

2.0 

1.0 

4-49 

7(s) 

1.2 

3.0 

0.1 

4-50 

7(t) 

1.2 

3.0 

0.5 

4-51 

7(u) 

1.2 

3.0 

1.0 

4-52 

7(v) 

2.0 

3.0 

0.1 

4-53 

7(w) 

2.0 

3.0 

0.5 

4-54 

7(x) 

2.0 

3.0 

1.0 

4-55 

7(y) 

3.0 

3.0 

0.1 

4-56 

7(z) 

3.0 

3.0 

0.5 

4-57 

7(aa) 

3.0 

3.0 

1.0 

4-58 


4-4 

GENERAL DYNAMICS CONVAIR DIVISION 



GENERAL DYNAMICS CONVAIR DIVISION 


t 

VI 


T BAR X / T s 1 .200X10*®° L/R = 1.000X10*°° T BAR » / T * I.OOOXIO* 01 


.01000 



GDC- DDG- 67-006 



GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s I.BQOXIO* 00 L/R = I.OOQXIO* 00 T BAR » / T = 5.000X10 



GDC-DDG-6 7-006 



general dynamics convair division 


T BAR X / T « t .800X10*°® l/R * 1.00QX10* 00 T BAR r t T * i. 000X10 



GDC- DDG-6 7-006 


GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X ✓ T = R.GOQXIQ * 00 L/R * I.OOOXIO * 00 T BAR Y / t a I. QOOXlO 



GDC-DDG-67-006 



GENERAL DYNAMICS CONVAIR DIVISION 


* .OOTOO 


c 

L 

A 

J .00900 

s 

e 



GDC-DDG-67-006 




4-10 

GENERAL DYNAMICS CONVAIR DIVISION 


t bar x / t * r.oooxiO * 00 l/r » j.oooxio * 00 t bar y / t * t.oooxio * 00 



GDC-DDG-67-006 


4-11 

general dynamics convair division 


T 0*M X / T I 1.00QX10* 00 l/« « 1.000XJ0* 00 T BAR » / T « 1.000X10 01 


* 

hhhhh 


p 

||g| 

Si 

”“T 

— 

— 

— 

■i 

■■1 

■HI 

PI 

Mil 

PI 

Ml 

■i 

11 


■■■■■I 

I'Hii i nm.M 

p 

p 

pi 

■ I 







■PPI 

PPPI 

PI 

■1 

■i 

■I 


■■■■■■HI 

p 

|PPP 

jggjjp 

ppi 

1 







1 







hhhhi 

P 

P 

p 

pm 

P| 













•00900 

■mmhh 

p 

p 

p 

Mi! 





hmbmh i 

mmm\ 

PI 

PI 

HI 






mmhhm 

p 

p 

p 


Mil 


_ 

"" 










— 


■ 

PPPPP 

P 

p 

up i 















PPPPP 

p 

pi 

pi 














*00000 

s 

l 

0 

* .00700 

B 

A 

R 

/ .00000 
e 

L 

A 

* .00000 

1 

c 

c 

.00400 

•00000 

.00000 

•00100 

•0 

rwnMHH 

PPPPPP 


ppp 

jpp: 













— 


HHP 

p 

PPP 

PI 

Mi 




pppppppppp; 

mmmm\ 

HM 

PPM! 

PPPI 

PPI 


_ 

pi 


wfiirnfn 

l« 








HI 

ppi 

PPI 

PPI 





P 

pppp 

1W 
















P 

pppp 

vm 










zd 





■HMHi 

p p 

pppp 

ppppp 





pppppppppp 

PPPPPPP' 

— I 

PPPP! 

PPPI 

g mm 





P 

ppn ^ 

nr 

ipp 





ppppppppppp 

PPPPPP 

■■■! 

PPPI 

ppi 

_J 





PPPPP 

pppp 

~T" 

\ 





pppppppppp 

PPPPPP 

P 

PP 

ppi 

PP 




ppppppppppp 

■■■■ 

pppp 


\ 





HHHMHMMM 

PPPPPP' 


PPM 

pp 







ppppp 

ppi 

\ 





ppppppppppp 

PPPPPP 

PPPPP 

PPM: 

ppi 

PP: 




pp p 

■ 

ppppjf 

pppM 






pppppppppp 

piPp 

PPPP 

PP 

■HBI 

PHi 





ppppppppppp 

■ ■ 

Hi 

ppp 

1 















imi 


PPP 

p 









ZZ 






p 


i PPPPPP 

ippp 

ip 

L. _ 




pppppppppip 

P PP 

PPPP 

PPPP 

PPP 

HP 




SI 


wmmm 

[pppp 

[ppp 


ilP 








: 


PP' 

9 

9 



IMP 

ippp 

mm 

n 










PP 


p 

\mmmum 

i S 

IPM 

mm 

X 




ppppppppppp 

PPHPMP 

PPPPP 

PPP 

PP 

Pi 




pp| 

IHMBI 

■ 



~t 


!■ 

■ 










p 

IHHB1 


< 

>^1 

“V 

V 

5 \ _ 












p 

\ummm 

PM 

I *4 

b 



PKilliSiPMHMMf 

1 PPPPPP 

PIMP 

PPM 

PP 





PR 

\mmmm 

ipppp 


*3 


lE 



ppppppppppp 

PPPPPP 

p P 

PPM 

PPP 

PR 







i „ 



mi 

:■ 

■ 

PPPPIPPPPPP 

PPPPPP 

IPPPP 

PPP 

HP 





pppppppp 

PPPPPP 

ppppp 

PPP 


'■ 

IBM 

r 

ppppppppppp 

IHHHHBP 

PPPP 

PPP 

ppp 





ppppppp 

PPPPP 

pppp 

PPP 



1 

m 

m 

r 

ppppppppppp 

PPPPPP 

PPPPP 

PPPI 

PP 

RPR 




pppppppp 

PPPP 

ppp 

PPP 

PPP 

IPP 

mi 

in 

p 

pppppppppp 

PPPPPP 

PPPP 

PPP® 

P 

RPR 




pmwTiiwim 

PPPPP 

pppp 

PPP 

IMI 

El 

e: 

>.w 

■ 

ppppppppppp 

PPPPP 

PPPP 

PPP 

PP 






PPPPP 

pppp 

PP 

pp 

IM 

!■ 

INI 

■ 

ppppppppppp 

PPPPP 

PPPP 

IPPP 

PPP 





ppppp 

PPPPP 

pppp 

PPP 

pp 

jpp 

!■ 


M 

pppppppppp 

PPPPPPP 

PPPP 

IMP 

ppp 






pppppppp 

PPPPPP 

ppp 

PPP 

ppp 

ipp 



P 

pppppppppp 

PPPPPPP 

pppp 

IPP 

PPP 

IPH 

IPP 




PPP 


IBM 

pp 

jpp 



■ 

IZHHHH 

PPPPPPP 

PPP 

IPPP 

PP 

L_ 





T 

r~~. 

i 

EI 

EZ 



e: 

Pks^PHPRPM 

PPPPPP 

IPPP 

IPPP 

PP 

IPH 




r 


PPPPPPPP 

ppppp 

IPP 

PP 

IPP 




-Figure 6(g) - Critical Axial Compress! 

for General Instability 
‘ Cylinders With Eccentric 

! Orthotropic Stiffening 

STRINGERS OUTSIDE 

RINGS INSIDE 

>Ov 





. .h 




o n v^.. 









in 




























VI 




















IJ S 











Has 















1 - - - 1 1 - L ■ 1 I I I I I 

10.0 100.0 1000.0 

L / RHO X 


GDC-DDG-6 7-006 


4-12 

GENERAL DYNAMICS CONVAIR DIVISION 


ft 

I 

G 

M 

A 

B 

A 

R 

/ 

C 

L 

A 

ft 

T 

I 

C 

c 


•01000 


.00000 


•00000 


•00700 


.00000 


•00500 


•00400 


•00500 


•00000 


•00100 


•0 



§ 

0 

1 

S 

? 

o 

o 

o> 





4-13 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T « 3.000X10 400 i./R a I.000X10* 00 T BAR Y / T * I.OOOXIO* 00 



GDC-DDG-67-006 


4-14 

GENERAL DYNAMICS CONVAIR DIVISION 


•OftOOO 


T BAR X / T = ! *200X10 


L/R s 2.000X10 


T OAR T / T s 1.000X10 


♦00 


♦00 


-01 



GDC-DDG-67-006 



4-15 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s |.800X»0 +00 - l/R * B.OOOXIO 400 T BAR T / T s S.QOOXIO' 01 

.otooo 



900 “i 9 -Daa-D(E) 


4-16 

GENERAL DYNAMICS CONVAIR DIVISION 


s 

I 

Q 

M 

A 

B 

A 

R 

/ 

£ 

L 

A 

S 

T 

I 

C 

E 


.01600 


•00000 


•00000 


•00700 


•00000 


•00500 


.00400 


•00300 


•OOtOO 


•00100 


•0 



GDC- DDG- 67- 006 








4-17 

DYNAMICS CONVAIR DIVISION 



T OAR * / T s 8.00QX10* 00 U/R * 8.000XI0 400 T BAR T / T « l.OOOXlfl’ 01 


•OIOOO 



GDC- DDG-6 7-006 




4-18 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T a B.OOOXIO * 00 t/R a B.OOOXIO * 00 T BAR Y I T a *. OOOXtO 

.01000 


.00000 

.00000 

•00700 

.00000 

.00500 

•00400 

•00900 

.00000 

•00100 

•0 

— 















































































m 

■ 

















m 

!■ 












































IB 

a 










HHHHHH 

IHHI 

!^H 

H 

hh 



■ 

HHHHHHH 

BiHBI^I 


■ 


hr 




H HH 

IHH 

■H 

H 

H 

HI 









_2 

u 

m 

IS 

MHUH 

1 
















U 













“ L 


Ll 

is 

m 

. \ \A -°v 



. 












11 

n 

■ 















m 

w 

■ 

\ V\n * 












LL 




s 

V \ 

& . 
















M 
















■1 

m 

■niAnHa 






m 

— 









□ 







3 









n 

hhmBII 

















■H^H 






















■ 

m 

-• 










hhhhhimh 





3 

3 










HHH^HH&\v\HH 








■ 

■■■1 

■HI3 

■Hi 

HH 

Hi 




HHHHHMvVH 






~ 


HHHHHHH 

■HHHI 

■HH 

HH 

H 





HHHHHgHHg^^ 




1 

m 

■ 

■ 






i 




hhhhhSh^m 




9 

9 

3 

a 










HHHHHuW^ 

V 



m 

M 

■ 

a 

RWi i Wi 10 1 H 

■ 

■ 

HI 

HI 

Hi 




HHHHnw 

V \H 



9 

m 

3 


■ ■ 

■■1 

■Hi 

HH 





hhhhhhhhkw 

\\\H| 






■ 


■ ■ 

hhh 

Hi 

HI 






H^HHfil 

ISKKB 








■ ■ 

■ 

■HI 

HH 








AIKS 







■ HH 

■HHH 

HHHI 

■Hi 

HH 

■H 




Ml 

HHK! 








HHHHHHH 

HHHH 

HHH 

■H 

HH 

_ 











- 

BHBBHHBBH | 

HHHH 

HHH 

HHI 

■H 

r 






hlsr 






Fi mire 6(n) - Critical Axial Comnre 

as 

.1 AfU 

HH 

BBBBS 







° 

for General Instability 
Cylinders With Eccentri 
Orthotropic Stiffening 

STRINGERS OUTSIDE 

RINGS INSIDE 

^ n 

HHHH 

HHH 






— 

m 

H ■ 

HHH 

hh : 

\^i 





c 

HH 

HH 

HH 

ss^ 

r 






HHH 

■MS 

HH 

hh 













H 


— 





‘ 

■ 


■ 

! 















H 

■ 











HHH! 

H 

HI 

HH; 


■ 

■; 

a 




PH 

HI! 









- 


Hi' 

■ 

i 

■■■■■I; 

■■■ 


■ 

■ 











■ 

3 

■ ■ 

Hi 

■ 


■ 













a 

3H HH 

i HI 

h 

■ 

H 





. 







Hi' 


HHHHHH 

■HI 

HH 

H 

■HI 

Hi 












□ 

n 


ioo.o iooo.o 

t- / RHO X 


GDC-DDG-6 7-006 



4-19 

GENERAL DYNAMICS CONVAIR DIVISION 


ft 

1 

0 

M 

A 

8 

A 

R 

✓ 

E 

L 

A 

ft 

T 

1 

C 

C 


T BAR X / T * 2*000X10*®° t/R = 2*000X10*°° T 8AR Y ✓ T * I.OOOXIO* 00 



10*0 


1000*0 


GDC-DDG-6 7-006 



GENERAL DYNAMICS CONVAIR DIVISION 


.01000 


T BAR X / T : 5.000X10 


L/R s 2,000X10 


T BAR T / T s 1.000X10 


I 

to 


O 


♦00 


-01 



GDC-DDG-6 7-006 



4-21 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 





4-22 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T " 3.000X10 


♦00 


L/R s 2*000X10 


♦00 


T OAR Y / T s 1*000X10 


♦ 00 


•01000 


$ 

I 

0 

H 

A 

B 

A 

R 


E 

L 

A 

S 

T 

1 

C 


•00000 


•00000 


•00700 


•00000 


•00300 


•00400 


•OOS0O 


•00200 


•00100 






r, ‘ 

■** 

'* 


““ 








MM- 









































































































— 

























































_ 









p 


\ '\ \ 














Zj 

z 

T 

t 

T \ 

\ 















Jt 

A 


L. 1 

















V \ \ 














z. 



L 


Y " 1 
















Y 

v \ \ 

V to 


*“ " 













3 

\ \ \ 

^ O 

















V V \ \ 













. 




V \ \ ^ V ~ 


































\ \ 'Ox \ \ 











... 






\ \ J o \ \ 

















\ | \ V \ 

















\ \ \ \ 

















V °o \ Y 


















\ 
















V \\\ 

















. VY 

Y\ 












1 













1 



i >< | 


l 1 






YVT 







Figure 

6(r) - Critical Axial Compi 
for General Instabi 

ression 

\ 

\v. 







lity in 

A 

vX\\ 

\ 







ntric s *“”" 


\\\ 

v\ 






i uyjLinaers nun wwc 



\N 

sSY 







Orthotropic stiffen 
STRINGERS OUTSIDE 

mg 


. \ 


W 



























, 





KimiS XIMaXLl** 














— 

— 

— 

— 

— 













■ r 



~ 

— — 

— — 





— 

— : 



















-r ~ 1 _ 



■ ■ 

1111™ 

WsM 

Pi 

2 

~ — 

— 

— 

— 

■ 

■ 

■ 

— 

-- 


~ 


















— 

— 

— 

— 















’ - Lr - 

■* 


— 

n n 




— 

~ 

— 

— 

— 







— 

— 


a 


I- / RHO X 


1000*0 


GDC- DDG- 6 7-006 



4-23 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s 1.200X10 


♦ 00 


L/R s 3.000X10 


♦00 


T BAR i / T * 1. 000X10 


’01 



GDC-DDG-67-006 



4-24 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T = 1.200X10* 00 L/R * S.OOOXtQ* 00 T BAR V / T s S. OOOXIO 



GDC-DDG-67-006 



4-25 

GENERAL DYNAMICS CONVAIR DIVISION 


s 

l 

0 

M 

A 

B 

A 

R 

/ 

E 

L 

A 

3 

T 

I 

C 

C 


t BAR X / T s 1.200X10*°° L/R * B.OOQXtO* 00 T BAR t / T a i.OOOXlO* 00 

>01000 ■ __ 




■■■bii 

■HI 

■■i 

Mil 

“i 















ismuhi 

■HI 

■i 

■II 

Ml 













.00*00 

.00*00 

•00700 

•00*00 

•00900 

•00400 

•00300 






ZJ 


















31 













■hhhh 

— i 

HI 


HI 

■II 













hhhh 

■^^■i 

HI 

mm 

■■1 

■ 




■^■■■■HB 

MMM 

■M 

mh 

■ 

■ 













MBBMHIIM 

mmmm 

■M 

■ 















■■ ■ 

B ■ 

■ 

■ 

BM 















■■■■ 

gglM 


BM 





■■■■■■■1 


M^| 

■■■ 

■ 

■HI 





■ ■ 

■MB 

■Mi 

Mi 





■MI^^^IBBI 

■ 


■M 

H 

■n 




V 


V \ 

zsz 

■ 

i 













\ 

v 

xz'srr \Zj 









33 






\ 

\ ~o \ 

<»> 

r. 





!■■■■ 

■■M 

^M 

■■ 

■i 




V 


\ A J 

L* 





IHMMMMM 

■ ■ 


mb 

MB 






\ 

<<v \ 
















MMSmB 

L^^BBB 







33 






V . 

r^T 

m 





■■■■■■Mi 

■ 

■■ 

■M 

■■ 

■i 







■mi 

ZL 














■ ■ 

M^BM 

- 

31 

33 














■ ■ 

Mi 

■aw 















■MBBMHM 

mbbbb 

IMW 

3“ 

S3 


■ 

■ 











■■1 

HMI 

M 

mm 













MHBBHHIB 

MBBBB 

■EOM 

IM 



■i 

■ 










■Mpiyffj r jgSMB[ 

■ ■ 

mm 


MB 




■ 

■■■i 

■MM 

■■■ 

■■1 

■■ 

■ 




■■■■■■■ 


mbm 

331 


ra 

kMI 



■ 

h 

■MB 

■ 

MB 





■■■■■■Mi 

■ ■ 

■ 

MK 1 

BM 

NS 

M 



■ 

.■■Hi 


■me 

■■ 





■ 

■ ■ 

■ ■ 

BBBi 

W 

BK 




■ 

'■HH 


wmm>. 

;■■ 





■ ■ 

■ 

■MB 

■Ml 

MSI 

W 

K1 



■ 

IHH1 


■M 

■■ 

mm 




■■ ■ 

■ ■ 

MBM 

BBS 




■i 


■ MB 

■■■ 

■M 

MH 

IHB 





■ 

■■■■■ 

■ 

■M 

Ml 

Kt 

2W 

■ 

■ 


I 

/ \ r* 


<i A 

i 

r^“ 

■« r 

. 



1 j 

■MMM 

■M 

■Mi 

BM 

hk 

IS 

\i 

M 

Figure 6 

vu; - urilicaj. axial v.omprea 
for General Instabilit 
Cylinders With Eccentr 
Orthotropic Stiffening 

STRINGERS OUTSIDE 

•sion 

■ 

■mb 

MK 

■M 

rz 

n 

v 

A 


y In 

■Ml 

■ 

■ 

I^M 



V 

K 


k l ft 




U- 



K 

Kj 











•ootao 

•00100 

•• 



■1 

■■■ 

mm 





j 



iflBBBBB 

■ 

MM 

n 

■ 















R 



T NaS TNSTDfc 





























■ 

■ BM 


IH 

■ 

■ 




■ 

■ 

■ 

■M 

mh 

■ 




■ 

■■ 

H 

IH 

■ 















BBIM 

■ 

■H 

IMH 

■ 

■ 













MBB 

■ 

MH 

!■■ 

BH 














MHB^^^MMI 

MMMM 

!■■■ 


1 





■ 

■ 

!■■■ 

■ 

IH 

■ 




BBBMBMBBM 

■■Mi 

■M 

■ 

■ 














■ 

MM 

■Mi 


M£gj 














■BBBI 

■M1M 

1— 

■ 

■ 














MMMMM 

■M 

IMHB 

■ 

■ 

■ 




IBIilHHBH 

IHHI 

IH 

IIBH 

IMI 

;M 




BBBMMBHM 

■■■Mi 

■M 

■M 

BB 

IBB 


amm 


IHHHH 

\mmmmm 

IH 

IH 

IBM 

IBB 





100.0 

l* / RHO X 


so.o 


ft 000.0 


900-i9-Daa~0CE) 



4-26 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T S 2.000X10*°** L/R s 3.000X10* 00 T BAR Y / T = 1.000X10 



GDC-DDG-6 7-006 




4-27 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T » 2.000XIQ* 00 L/R = 3.000X10* 00 t BAR T / T = S. OOOXIO 



GDC- DDG- 6 7-006 



4-28 

GENERAL DYNAMICS CONVAIR DIVISION 


1 

0 

M 

A 

B 

A 

R 

/ 

E 

L 

A 

S 

T 

I 

c 

E 




9oo-i9-Daa-oaD 







4—29 

GENERAL DYNAMICS CONVAIR DIVISION 


s 

1 

6 

M 

A 

B 

A 

R 

/ 

E 

L 

A 

S 

T 

I 

C 

E 



GDC-DDG-6 7-006 




4-30 

GENERAL DYNAMICS CONVAIR DIVISION 


s 

X 

G 

H 

A 

B 

A 

R 

/ 

£ 

t 

A 

ft 

T 

1 

C 

£ 



GDC- DDG- 6 7- 006 


4-31 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / t * S.000X1Q* 00 * S.000XJ0* 00 T »*« t / T « I.OOOXIO* 00 



GDC-DDG-6 7-006 



4-32 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s 1.200X10* 00 L/R * I.OOOXlQ* 00 T BAR t / T * i.OOOXlO* 01 

. 0*000 



GDC- DDG- 6 7-006 



4-33 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s I.ROOXIO * 00 L/R * I.OOOXIO * 00 T BAR Y / T « J.OOOXIO * 01 



GDC-DDG-67-006 



4—34 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s 1.20QX10* 00 L/R s 1.000X10* 00 t bar T / T S 1.000X10* 00 

.01000 



GDC- DDG- 6 7- 006 





4-36 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s B.OOOXIO* 00 U/R * 1.000X10* 00 T BAR T / T * 5.000X10 

.01000 



GDC-DDG-67-006 




4-37 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T = 8.000XI0* 00 1/H = ».000X»0* DO T BAR T / T * I.OOOXIO* 00 



GDC-DDG-6 7-006 


4-38 

GENERAL DYNAMICS CONVAIR DIVISION 


s 

1 

G 

M 

A 

8 

A 

R 

/ 

E 

L 

A 

S 

T 

1 

C 

E 


•OlOOO 


•00900 


•00900 


•ooroo 


.00000 


•00500 


•00400 


•00900 


•ootoo 


•ootoo 


•0 



fl 000.0 


GDC-DDG- 6 7- 006 




4-39 

' GENERAL DYNAMICS CONVAIR DIVISION 




4-40 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T * J.OOOXIO * 00 , ' L/R - i.OOOXlO * 00 T BAR r / T = I.000X10* 00 



GDC- DDG- 67-006 


4-41 

L DYNAMICS CONVAIR DIVISION 


T IU I / r < t.tOOXlO * 00 


U/« » *.000Xlfl t00 


T •*« t t T • l.OOOXIO 


-01 


.01000 

SEE 

MSI 

wamm\ 

SHI 

s 

5! 

s; 

3 

_ 

— 

— 




— 

— 

— 

— 


u 

8 

•ootoo 

SEES 

sssi 

ssn 

SHI 

■SSI 

MSI 

= 

BS! 

BBI 

wm\ 

MHMI 

S! 

3! 

MM] 


— 

— 

— 




— 

— 

— 

— 


B 

■ 

1 

•00000 

3353 

p 

■ 

pm 

■ 

bll 

BBi 

■ 

— 

■S 

■l 

— 




— 

— 

— 

— « 

— 

— 

s 

1 

0 

JJ « OOTOO 

0 

A 

R 

, ,00000 
e 

L 

A 

* .00500 

t 

c 

t 

•00400 

•00300 

•ootoo 

•ootoo 

•i 




■Ml 

CIV 

mm 

— 

Mi 

— 

kl 




r 

■ 

■■ 



— 




■11 



m 


T 




B 

BB 




, 

p 

■ P 

p 

All! 

l^P 

u 

' 

: 

□ 

BBSs 

bb 

BBB 

BBI 

BH 








■IV 




■i 

■.Us 

BBI 

BBB 

MB, 

p 




. 




Kill 

AM 




□ 

IUS 

■ 

BBB 

MB; 

BH 




-J . . 

r — 




■111 

■HU 

M 

MSI 


3 

3 

\° 


■ 

■ 

Sj 

p; 

m 

— 

— 

— 




Bill 

■\u 

■nil 

1BH 

mm 

ASS 

BB 

ESI 

m 

i= 


pip 

p 

■ 

«■ 


— 

— 

— 



ppp 

BkVt 

kml 

1W 

H 

W1 

m 


p 

isuaiiao 

■■■■■ 

p 

pp 

■ 

P 

"" T 



p p 

p p 

p p 

pv 

IVI 

s 

R 

m 

» 

■M3PH 

BB 

BBB 

BB 

BB 

f-J 





p 

p p 

p p 

ppj 

111 

n 

|g 

S 

P 

MLWHPH 

PPBP 

BBB 

BB 

BB 

— 



m 


■ 

p 

p 

111 


a 

s 

IP 

H3HBBBBB 

PPP 

BBB 

BB 

BB 

■B 

^ ( 


mmmr 





mm 


m 

'31 

PMIWHPP 

BBI 

BBB 

BB 

BB 

__ 








in 

VI 



k 

ppjpp 

bbbb 

BiltTTB 

BB 

^B 

— 









IV 

UK 

LZ 

V 

r 


B 

BBB 

BBB 

BB 










■% 


m 


\ > 

ppppppiPi 

BBB 

BB 

B 

BB 

BB 











iUK 



Cf#pppp|^HKB 

BBBB 

BBB 

BBI 

HB 




_ 






11 

IM 

im 

■ 

nmw 

BBBB 

BBB 

BB 





. 






IS 

■i 

LZ 

IS 

snmv 

I33333 

BBB 

BBI 


_ 











IKI 

IM 

Iftfi 

m 

1 ^BBB2!HK9H^ 

^ 



_ 













IBI 

:v 

r_ 

K 


L%KBBB 










ppp 

■ 

|PP 

IS 

\a 

laa 


II TMlilV 

BBi 







IB 


p 

ppp 


IPS 


> 

!«• 

£ 

iKJn 








a 

fSi l il iiiii'iTFifg^maO^ 

p 

ppp 

T~ 




IK1 

IK 

ismkSSIk? 

IKK^I 

iBBi 




_ 





s 


p 

pilpyy 

!■■ 






ST 


!BB 






IE 


p 

p 








IlkS^ 









s 

p 

!B 







!BsK! 

I5SSM 








— 

— 


— 






- 5sn -» 

— it, ^ 

i:sa*.s$s 

: ?BB 






Figure 7(j) - Critical Axial t-ompressi 
for General Instability 
‘ Cylinders With Eccentric 

Orthotropic Stiffening 

I STRINGERS INSIDE 

I RINGS INSIDE 

on ^*= 









in 



IjjJ^ ^P 





■PPPP 

BBBB 

BBB 






■ppp 

bbb 

BBB 

pp 

■^5 




pp ppp 

bbbb 

BBB 

BB 

BB 

llife 



JBBBBBB 

bbbb 

BBB 

BB 

BB 

EZ 



5 

Jbpbppb 

bbbb 

BB 

BB 

BB 

tz 



IE 

pp 

BB3 

IBBB 

BIB 

IPi 

IIP 



i9 

■bb 

BBBB 

IBBB 

BB 

BH 

1 



ieI 


I ■ * * . ■ ...... 

, 0 . 0 100 . 0 1000. 0 

L / RHO X 


GDC-DDG-67-006 



4-42 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X ✓ T = 1.200X10 


♦ 00 


L/R = 2.000X10 


♦ 00 


T BAR Y t T = 5.000X10 


-Cl 



GDC-DDG-6 7-006 



GENERAL DYNAMICS CONVAIR DIVISION 



GDC-DDG-67-0O6 



4-44 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T = 2,000X10 


♦00 


L/R = 2*000X10 


♦00 


T BAR V / T s 1,000X10 


-01 


5 

1 

0 

M 

A 

B 

A 

R 


C 

L 

A 

S 

T 

I 

C 


•01000 


•ootoo 


•00000 


•00700 


•00000 


•00900 


.00400 


.00900 


•00200 


.00100 


«0 



10.0 


100.0 

L ✓ RHO X 


1000.0 


GDC-DDG-67-006 


DYNAMICS CONVAIR DIVISION 


I 

& 


T BAR X / T s R.OOOXIO * 00 UR * ttOOOXtO * 00 T BAR t / T » J.OOOXIO" 0 * 



GDC- DDG-6 7-006 



4-46 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T - 2.000X10 


L/R = 2.000X10 


T BAR Y / T s 1.000X10 


♦ 00 


♦00 


♦ 00 



GDC-DDG-67-006 



4-47 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T * S.000XI0* 00 L/X » t.OOOXtO* 00 t •** t / T ■ I. OOOXIO 0> 



GDC— DDG-67-006 



4-48 

GENERAL DYNAMICS CONVAIR DIVISION 


•01000 


T BAR X / T s 3*000X10 


L/R = 2.000X10 


♦00 


♦00 


T BAR Y f T 


3*000X10 


-01 


t t i 1 1 — — r — rh i r 

Figure 7(q) - Critical Axial Compression 
for General Instability in 
Cylinders With Eccentric 
Orthotropic Stiffening 


STRINGERS INSIDE 



•0 


10*0 


100.0 


1000.0 


GDC-DDG- 6 7-006 



4-49 

GENERAL DYNAMICS CONVAIR DIVISION 





4-50 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T * I.BOOXIO * 00 UR = J.OOOXIO * 00 T BAR » / T s J.QOQXlo ’ 01 

. 0*000 



GDC- DDG-6 7-006 




4-51 

GENERAL DYNAMICS CONVAIR DIVISION 



9oo-£9-£xia-oaD 



4-52 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T = I.800X10* 00 L/R = 3.00QXJ0* 00 T BAR Y / T * l.OOOXlO* 00 


.01000 



GDC- DDG-6 7-006 








4-53 

general dynamics convair division 



GDC- DDG-6 7-006 



4-54 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T s 2.000X10 


L/R s 5. 000X10 


T BAR T / T s 5.000X10 


Klim 


iwma 





■■ 


■■ 














■Ufil 





\m*2m 

mmm 


i»i 


iw 


Figure 7(w) - 


Critical Axial Compression 
for General Instability in 
Cylinders With Eccentric 
Orthotropic Stiffening 

STRINGERS INSIDE 


RINGS INSIDE 


mil 


BMgm msiiawgqll 


(KsSisiiii 


100.0 

t / RHO X 


GDC-DDG-67-006 



4-55 

GENERAL DYNAMICS CONVAIR DIVISION 


T IAR X / T * t, 000X10* 00 L/R * R.OOOXIO* 00 T BAR Y / T ■ 1,000X10 



GDC— DDG-67-006 



4-56 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T = 3>000X10*°° L/R = 3.000X10*°° T BAR Y / T = 1,000X10 



GDC- DDG- 67-006 


4-57 

GENERAL DYNAMICS CONVAIR DIVISION 


t BAR X / t * 3.000X10 


♦ 00 


L/R * 9*000X10 


♦00 


T OAR T / T ■ 9. 000X10 


*01 


.01000 


9 

1 

0 

M 

A 

B 

A 

R 



100.0 

t / RHO X 


1000*0 


GDC-DDG-67-006 


4-58 

GENERAL DYNAMICS CONVAIR DIVISION 


T BAR X / T » S.QQQXiO* 00 t/R = 3.000XI0* 00 T BAR V / T s I.OOOXIO * 00 


•OIOOO 



GDC- DDG-6 7-006 



GDC-DDG-6 7-006 


SECTION 5 
ELASTIC CONSTANTS 

The digital computer program of SECTION 6 includes an option which 
allows for the input of elastic constants and eccentricity coupling con- 
stants. This feature was incorporated to provide the engineer with a more 
flexible analysis method than is given by the curves of SECTION 4. However, 
in order to make use of this capability, one must first compute appropriate 

values for the A. .'s, D. .'s, and C. ,'s. Some recommended formulas for these 
ij 13 13 

constants are listed in Table III. The tabulated formulas are simplified 
expressions suitable for application to most practical stiffened configura- 
tions (not including waffle-type structures). To be rigorous, more com- 
plicated expressions would be required. All of the given formulas apply 
only where the behavior is elastic. For cases where the buckling stress ex- 
ceeds the proportional limit of the stress-strain curve, the stringer (or 
corrugation) and skin stiffnesses should be modified to reflect tangent 
modulus influences. 

In addition to the elastic constants and eccentricity coupling constants. 

Table III also lists recommended formulas for the computation of the quantities 

t , t , p , and p which must be computed if one is to make use of the curves 
x y x y 

given in SECTION 4 or the RATIO option of digital computer program 4267 (see 
SECTION 6). 

To fully understand Table III, it is helpful to note that the A *s 
and D„'s arise out of mathematical integrations involving the distribution 
of the composite wall material about the appropriate centroidal surface. 

Note that the centroidal surface has a curvature of its own. Therefore, 
the related material distribution is equivalent to that which exists about 
the centroidal plane of the flat plate obtained by unfolding the composite 
circular shell wall into a flat configuration. All influences of curvature 
in this regard are inherent in the basic shell equations into which the 
A. .'s and D. . »s are substituted. 


5-1 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 


Table III applies only to cases where no buckling of the isotropic 
skin panels and no local buckling of the stringers occur prior to over- 
all instability. In addition, it is assumed that the stringers are spaced 
sufficiently close together to justify the assumption that all of the skin 
material is fully effective. In cases where these several conditions do 
not prevail, it becomes necessary to introduce effective-width concepts to 
modify the information given in the table and the notes which follow it. 
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Notes for TABLE III 


(a) For convenience, all of the non-corrug.ited configurations shown here only 

depict the case of external stringers. In addition, only internal rings 

are shown. However, all of the formulas for t , t , P , P , the A. .’s.and 

* x’ y’ x’ y’ lj ’ 

the D. .'s apply equally well for other orientations of the stiffeners (inside 
vs. outside). On the other hand, the signs of the and C 22 values depend 

on the stiffener locations (z^ is positive for internal stringers and negative 
for external stringers; z is positive for internal rings and negative for 

y 

external rings). 

(b) The quantity t x is the wall thickness for a monocoque circular cylinder of the 
same radius as the middle surface of the sti f f ened-cylinder basic skin, and of 
the same total cross-sectional area as the actual composite stiffened wall in- 
cluding both skin and stringers. The cross section referred to here is obtained 
by passing a plane through the cylinder, normal to the axis of revolution. 

(c) The quantity A is the cross-sectional area of a single stringer and does not 

s 

include any of the basic cylindrical skin. The cross section referred to here 
is obtained by passing a plane through the cylinder, normal to the axis of 
revolution. 

(d) The symbol (SdJ is used to denote the total peripheral length of the corruga- 
tion center-line for the wave-type cross section obtained by passing a plane 
through the cylinder, normal to the axis of revolution. Hence, (2d^)> 2itR and 
the total area for the stated corrugation cross section may be taken equal to 

(Sd.)(t ). 

1 c 
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Notes for TABLE III (Continued) 


(e) The quantity A is the cross-sectional area of a single ring augmented by 

r+s 

an effective width of skin which is considered to essentially behave as part 
of the ring. The user must apply engineering judgement in selecting the effective 
width to be used here. The cross section referred to lies in a radial plane. 

(f) The quantity A f is the cross-sectional area of a single ring and does not include 
any effective width of skin. The cross section referred to here lies in a radial 
plane. 

(g) The quantity a is the spacing between rings (assumed to be uniform throughout 
the structure). 

(h) The quantity T is the local longitudinal centroidal running moment of inertia 
for the flat configuration obtained by unfolding the entire Composite circular 
shell wall. For example, consider the case of a cylinder having a local wall 
cross section of the type 
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Notes for TABLE III (Continued) 



After computing 1^ ^ for the cross-hatched area shown, I is found as follows: 



(i) The quantity I is the centroidal moment of inertia for the cross section of 

r+s 

a single ring augmented by an effective width of skin which is considered to 
essentially behave as part of the ring. The user must apply engineering judgement 
in selecting the effective width to be used here. The cross section referred 
to lies in a radial plane. 
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Notes for TABLE III (Continued) 


(j) The quantity I r is the centroidal moment of inertia for the cross section 
of a single ring and does not include any effective width of skin. The 
cross section referred to here lies in a radial plane. 

(k) The term (l-v 2 ) has been omitted from the formulas for D^ and D 22 since 
the specified configurations provide incomplete restraint to the related 
anticlastic bending (see GLOSSARY, Volume I Cl]). 

(l) The several values which have been set equal to zero all involve the 
available mechanisms for interplay that are dependent upon Poisson' s-ratio 
influences. For usual stiffened constructions, such interplay is of a 
minimal nature since it occurs primarily in the localized regions where 
the stiffeners intersect. However, when the stringer and ring spacings 
are very small (as for waffle configurations) and/or the contributions 
from stiffener rigidities approach that of the skin (also as for waffle 
configurations), non-zero values should be used for A^g, D^g, C^g, and Cg^. 
However, no recommendations are made here as to practical formulations to 
be used for the computation of such values. 
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SECTION 6 

DIGITAL COMPUTER PROGRAM 

This section presents the essential features of General Dynamics 
Convair digital computer program numbered 4267. This program was 
developed for the analysis of general instability in axially compressed 
circular cylinders having eccentric stringers and rings. The solution 
is based upon the theoretical considerations presented in SECTION 2. 

The output can be obtained in the form of automatically plotted data or 
single-point solutions, as desired. All of the buckling curves presented 
in SECTION 4 were obtained by using the automatic plotting option of the 
program. The input format is shown in Figure 8. Symbols are listed in 
Table IV. A detailed, card-by-card description of the input follows 
below. Runs may be stacked. For further information concerning this 
program one may refer to the detailed derivations and basic logic pre- 
sented in reference 15. 

CARD TYPE 1: One card per run. 

Enter PROBLEM IDENTIFICATION anywhere in columns 1-60. 

Alphanumeric characters. 

CARD TYPE 2: One card per run. 

Enter INPUT OPTION (RATIO or STIFF) in columns 1-5. 

This option permits the user to choose between 
alternative formats for cards other than types 1, 2, 
and 3. 

Enter NO OF CASES as right adjusted integer in 
columns 6-10 (15). 

Enter POISSON’S RATIO (v) in columns 11-15 (F5) . 

Enter E/G (Young's modulus 4 Modulus of Elasticity 
in Shear) in columns 16-20 (F5). 
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Enter LOWER in columns 21-25 (F5) . Output 

is listed for critical ■ (=m* ) and also for 

m= (LOWER C ) x m*. The value (LOWER C ) = 0.99 
a a 

will usually be suitable. 

Enter UPPER in columns 26-30 (F5). Output is 

listed for critical m (=m*) and also for m = 

UPPER C ) x m*. The value (UPPER C ) = 1.01 will 
a a 

usually be suitable. 

Enter LOWER C D in columns 31-35 (F5). Output is 

P 

listed for critical p (=p*) and also for p = 

(LOWER C D ) x p*. The value (LOWER C_) = 0.99 will 
P P 

usually be suitable. 

Enter UPPER C Q in columns 36-40 (F5). Output is 

P 

listed for critical p (=p*) and also for p = 

(UPPER C D ) x p*. The value (UPPER C_) = 1.01 will 
P P 

usually be suitable. 

Enter NRRHOY (number of R/p ratios to be included 

y 

in plots and/or tables which result from automatic 
sequencing operations) as right adjusted integer 
in columns 41-45 (15). Will be left blank when 
only point solutions are to be obtained. 


Enter TJ OPTION as right adjusted integer in 
columns 46-50 (15) . Always insert the number 1 
here except where 

(a) it is desired to eliminate the p con- 

xy 

tribution in RATIO-option point solutions 


or 


(b) it is desired to eliminate the Tl con- 

P 

tribution in STIFF-option point solutions. 
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In the latter two instances leave columns 
46-50 blank. 

Enter OPTION as right adjusted integer in 
columns 51-55 (15). Always insert the number 1 
here except where it is desired to eliminate the 
C contribution in point solutions. (For other 

X M 

types of solutions the program always assumes 
= 0) . For those point solutions where one 
chooses to ignore the influence, columns 

51-55 are to be left blank. 

Enter OPTION as right adjusted integer in 

columns 56-60 (15). Always insert the number 1 
here except where it is desired to eliminate the 
contribution in point solutions. (For other 
types of solutions, the program always assumes = 0) . 

For those point solutions where one chooses to 
ignore the influence, columns 56-60 are to be 

left blank. 

Enter A^ OPTION as right adjusted integer in 
columns 61-65 (15). Always insert the number 1 
here except where it is desired to eliminate the 
A 12 con * r; *-b u ti on in any type of solution. In the 
latter case, columns 61-65 are to be left blank. 

Enter as a right adjusted integer in columns 
66-70 (15). This input constitutes the number of 
refinement cycles used to improve the accuracy of 
final computed values. Whenever (BETA FACTOR) ^ 1.05, 
the value Ng=5 should usually be satisfactory. Each 
single refinement cycle essentially cuts the final |3 
screening increment in half. 
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Enter DUMP OPTION as a right adjusted integer 
in columns 71-75 (15). Insert the number 1 
whenever supplementary diagnostic output data 
is to be printed out. Otherwise, leave blank. 

CARD TYPE 3: One card per run. 

Enter MAX ORDINATE FOR PLOTS in columns 1-10 
(E10.5). Any of the following values may be 
inserted here: 


.0005 

.001 

.006 

.0006 

.002 

.007 

.0007 

.003 

.008 

.0008 

.004 

.009 

.0009 

.005 

.010 


This entry is left blank when no plots are to be made. 

Enter SCREENING CUT-OFF in columns 11-20 (E10.5). 

This is a cut-off value used in the minimization 
process and must be set greater than the output 
critical strain. Whenever the computed critical 
strain exceeds the SCREENING CUT-OFF value, the 
printed results cannot be believed. In such cases, 
one should increase the input SCREENING CUT-OFF value 
and rerun the program. In general one will have little 
interest in critical strains above .010. Therefore, 
this value will usually be a suitable choice for the 
SCREENING CUT-OFF. 

Enter MIN NO CIRCUMF HALF-WAVES in columns 21-30 (E10.5). 
This is the minimum number of circumferential half-waves 
considered to be permissible for non- axi symmetric buckle 
patterns. A value of 2.0 should usually be inserted here* 
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CARD TYPE 4 (RATIO OPTION): 

There will be NRRHOY/8 (rounded to the higher whole 
number) cards per run. However, when NRRHOY = 0, 
no cards of this type are required. 

Enter RRHOY (R/P ) values, with a maximum of 8 to . 

y 

a card (8 E10.5). 

CARD TYPE 5 (RATIO OPTION): • 

There will be one of these cards for each case. 

Enter CASE NO as right adjusted integer in columns 
1-5 (15). 

Enter OUTPUT OPTION as right adjusted integer in 
columns 6-10 (15). 

1 = Tables with no plots 

2 = Tables plus plots 

3 = Plots with no tables 

4 = Point solution 

Enter the thickness ratio (t /t) in columns 11-20 (E10.5). 

x 

Enter the (Length/Radius) ratio L/R in columns 21-30 
(E10.5). 

Enter the thickness ratio (t ^/t) in columns 31-40 (E10.5). 

Enter the geometric factor C^ x in columns 41-50 (E10.5). 
This factor is defined by equation (2-26) and is not 
used for point solutions. Hence, when OUTPUT OPTION = 4, 
leave this entry blank. The value of Cp x = .33 will 
usually be reasonable for practical skin-stringer com- 
binations. 
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Enter the geometric factor in columns 51-60 
(E10.5). This factor is defined by equation (2-25) 
and is not used for point solutions. Hence , when 
OUTPUT OPTION = 4, leave this entry blank. When it 
is desired that the computer employ the equation 

<w ■ (J) 

insert the left adjusted word AUTOMATIC in columns 
51-60. 

Enter the geometric factor C in columns 61-70 (E10.5). 
This factor is defined by equation (2-30) and is not 
used for point solutions. Hence, when OUTPUT OPTION = 4, 
leave this entry blank. The value C = .40 will 

r 

usually be reasonable for practical skin-ring com- 
binations. 

Enter the geometric factor C in columns 71-80 (E10.5). 
This factor is defined by equation (2-29) and is not 
used for point solutions. Hence, when OUTPUT OPTION = 4, 
leave this entry blank. The value = .425 will 

usually be reasonable for practical skin-ring com- 
binations. 

CARD TYPE 6 (RATIO OPTION): 

There will be one of these cards for each case. 

Enter STRINGER LOCATION (INSIDE, OUTSIDE, SYMMETRIC) 
as left adjusted word in columns 1-10. 
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Enter RING LOCATION (INSIDE, OUTSIDE, SYMMETRIC) 
as left adjusted word in columns 11-20. 

Enter BETA FACTOR in columns 21-30 (E10.5). This 
is a stepping factor used in the minimization 
process. That is, screening is performed involving 
p values computed from 



( 6 - 2 ) 


The value (BETA FACTOR) = 1.02 should be suitable for 
most applications. 


CARD TYPE 7 (RATIO OPTION): 

This card only required for point solutions (OUTPUT 
OPTION = 4). 

Enter the longitudinal slenderness ratio L/p^ in 
columns 1-10 (E10.5). 

Enter the hoop slenderness ratio R/p^, in columns 
11-20 (E10.5). 

Enter the artificial slenderness ratio L/p in 

r xy 

columns 21-30 (E10.5). 

Enter the eccentricity-dependent ratio C^/R in 
columns 31-40 (E10.5). 

Enter the eccentricity-dependent ratio C^/R in 
columns 41-50 (E10.5). 

Enter the eccentricity-dependent ratio C^/R in 
columns 51-60 (E10.5). 

Enter the eccentricity-dependent ratio C^/R in 
columns 61-70 (E10.5). 
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CARD TYPE 4 (STIFF OPTION): 

There will be one of these cards for each case. 

Enter CASE NO as right adjusted integer in columns 

1-10 ( 110 ). 

Enter the elastic constant A^ in columns 11-20 (E10.5). 

Enter the elastic constant A ^ in columns 21-30 (E10.5). 

Enter the elastic constant A^ in columns 31-40 (E10.5). 

Enter the elastic constant A in columns 41-50 (E10. 5). 

oo 

Enter the eccentricity coupling constant C^ in columns 
51-60 (E10.5). 

Enter the eccentricity coupling constant C pp in columns 
61-70 (E10.5). 

Enter R (the radius to the middle surface of the basic 
cylindrical skin) in columns 71-80 (E10.5). 

CARD TYPE 5 (STIFF OPTION): 

There will be one of these cards for each case. 

Enter the elastic constant D^ in columns 11-20 (E10.5). 

Enter the elastic constant D in columns 21-30 (E10.5). 

Enter the elastic constant D^^ in columns 31-40 (E10.5). 

Enter the elastic constant D in columns 41-50 (E10.5). 

Enter the eccentricity coupling constant in 

columns 51-60 (E10.5). 

Enter the eccentricity coupling constant in 

columns 61-70 (E10.5). 

Enter the overall length L (this is NOT the ring 
spacing) in columns 71-80 (E10.5). 
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CARD TYPE 6 (STIFF OPTION): 

There will be one of these cards for each case. 


Enter BETA FACTOR in columns 11-20 (E10.5). This 


is a stepping factor used in the minimization process. 
That is, screening is performed involving P values 
computed from 



(6-3) 


The value (BETA FACTOR) = 1.02 should be suitable for 
most applications. 


A sample input coding form is shown in Figure 9. 

The program output consists of a listing and/or plots depending upon 
the options selected. A sample output listing for 

INPUT OPTION = RATIO 
OUTPUT OPTION = 1 

is shown in Figure 10. Typical plots are given in SECTION 4. A basic 
flow diagram for the program is presented as Figure 11 and a Fortran listing 
of the program is shown in Table V. 
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TABLE IV - Program 4267 Notation 


Program 

Notation 

Report 

Notation 

Description 

TITLE 

■ 

Problem identification. 

RAORST 

Input Option 

RATIO or STIFF. 

NCASES 

— 

Number of cases. 

POISR 

V 

Poisson's ratio. 

EDSIG 

(E/G) 

Ratio of Young's modulus 
to modulus of elasticity in 
shear. 

CALPHL 

Lower C 
a 

See description of input 
for CARD TYPE 2. 

CALPHU 

Upper C a 

See description of input 
for CARD TYPE 2. 

CBETAL 

Lower C 

P 

See description of input 
for CARD TYPE 2. 

CBETAU 

Upper Cp 

See description of input 
for CARD TYPE 2. 

NRRHOY 

NRRHOY 

The number of R/p ratios 

y 

to be included in^ plots 
and/or tables which result 
from automatic sequencing 
operations. 

IETAP 

T) Option 
P 

See description of input 
for CARD TYPE 2. 

IC12 

Option 

See description of input 
for CARD TYPE 2. 

IC21 

Option 

See description of input 
for CARD TYPE 2. 

IA12 

A^^ Option 

See description of input 
for CARD TYPE 2. 

NR 

n r 

See description of input 
for CARD TYPE 2. 

IDUMP 

Dump Option 

ft 

See description of input 
for CARD TYPE 2. 
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Program 

Notation 

ORDMAX 

SCREEN 

CHtfMIN 

RDRHOY 

NCASE 

ITYPEO 

THICKX 

FLGRAD 

THICKY 

CRHOX 

CHX 

CRHOY 

CHY 

STRING 


TABLE IV - Program 4267 Notation 
(Continued) 


Report 

Notation 


Max Ordinate 
For Plots 

Screening 

Cut-Off 

Min no. circumf. 
half-waves 


Output Option 


t /t 
x 

L/R 


t /t 

y 


Px 




Stringer 

Location 


Description 

Uppermost grid line for 
plots. 

See description of input 
for CARD TYPE 3. 

Cut-off value used in mini- 
mization process. See 
description of input for 
CARD TYPE 3. 

R/P ratios to be in plots 
or tables. 

Case number. 

1 -= Tables only 

2 = Tables plus plots 

3 = Plots only 

4 = Point solution 

Thickness ratio. 

Length-to-radius ratio. 

Thickness ratio. 

Geometric factor defined 
by equation (2-26). 

Geometric factor defined 
by equation (2-25). 

Geometric factor defined 
by equation (2-30) „ 

Geometric factor defined 
by equation (2-29). 

Stringer location (INSIDE, 
OUTSIDE, SYMMETRIC). 
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TABLE IV - Program 4267 Notation 
(Continued) 


Program 

Notation 

Report 

Notation 

RINGL 

Ring Location 

BETAF 

Beta Factor 

FLRHOX 

L/P x 

FLRHOY 

x/p y 

FLROXY 

L/P xy 

C11R 

C n /R 

C22R 

C 2/ E 

C12R 

c 12 /R 

C21R 

C 21 /R 

All 

4 11 

A22 

A 22 

A12 

A 12 

A33 

A 33 

Cll 

C 11 

C22 

C 22 

R 

R 

Dll 

°11 

D22 

°22 

D12 

D 12 


Description 

Ring location (INSIDE, 

OUTSIDE, SYMMETRIC). 

See description of input 
for CARDS TYPE 6. 

Longitudinal slenderness 
ratio. 

Circumferential slenderness 
ratio for point solution only. 

Artificial slenderness ratio. 

Eccentricity coupling ratio. 

Eccentricity coupling ratio. 

Eccentricity coupling ratio. 

Eccentricity coupling ratio. 

Elastic constant. 

Elastic constant. 

Elastic constant. 

Elastic constant. 

Eccentricity coupling constant 

Eccentricity coupling constant 

Radius. 

Elastic constant. 

Elastic constant. 

Elastic constant. 
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TABLE IV 

- Program 4267 

Notation 



(Continued) 


Program 


Report 

Description 

Notation 


Notation 

D33 

°33 


Elastic constant. 

C12 

C 12 


Eccentricity coupling 
constant. 

C21 

C 21 


Eccentricity coupling 
constant. 

FLENGT 

L 


Overall length of cylinder. 
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FORTRAN CODING AND DATA FORM 
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GENERAL INSTABILITY OF ECCENTRICALLY STIFFENED CYLINDERS 4267 


SAMPLE PROBLEM 

INPUT NUMBER POISSONS LOWER UPPER LOWER UPPER ETA SUB P CI2 C2I A12 

OPTION OF CASES RATIO E/G C SUB ALPHA C SUB ALPHA C SUB BETA C SUB BETA NRRHOY OPTION OPTION OPTION OPTION 


RATIO I 0,300 2.600 0 

MAX ORDINATE 
FOR PLOTS 

- 0 . 

2.00000E 02 

CASE OUTPUT 

NUMBER OPTION T BAR X/T 

1 l 3.C0000E 00 

STRINGER LOCATION 

OUTSIDE 


L/RHO X 


CRITICAL 
LOWER CASE M 

CRITICAL 
LOWER CASE N 

3. QOOOGE 

01 

1.000006 

00 

5.68383E 

00 

4.250006 

01 

l.OOOOOE 

00 

5 • 72265E 

00 

6.00C00E 

01 

l.OOOOOE 

00 

5.66978E 

00 

8.60000E 

01 

l.OOOOOE 

00 

5.59234E 

00 

l . 21500fc 

02 

l.OOOOOE 

00 

5.51530E 

00 

1 • 75C00E 

02 

l.OOOOOE 

00 

5.45213E 

00 

2 • 45 000 E 

02 

1.53299E 

00 

6.37875E 

00 

3 • 50C00E 

02 

2.77085E 

00 

7.63813E 

00 

4 • 95000E 

02 

3.9570OE 

00 

8.29345E 

00 

7.0000CE 

02 

5.21559E 

00 

8.73586E 

00 

l.OOCOOE 

03 

6.645376 

00 

9.Q7263E 

CO 


990 I. CIO 0.990 1.010 

SCREENING MIN NO CIRC 

CUT-OFF HALF-WAVES 

1.00000E-02 2.00 

R/RHO Y RATIOS FOR AUTOMATIC SEQUENCING 


L/R 

T BAR Y/T 

C SUB RHO X 

•OCOOOE 00 

5.00000E-01 . 
RING LOCATION 
INSIDE 
R/RHO Y 
2.00000E 02 

3. 30000E-01 

CRITICAL 

LL 

LH 

SIGMA BAR /E 

SIGMA BAR/E 

SIGMA BAR/E 

2.599646-02 

2.56541E-02 

2.566086-02 

1.54430E-02 

1 • 5272 IE-02 

1.527766-02 

9.902496-03 

9.81826E-03 

9. 82342E-03 

6. 79592E-03 

6.75770E-03 

6.76227E-03 

5.15497E-03 

5. 139 17E-03 

5 • 143976-03 

4.17198E-03 

4.16858E-03 

4* 17285E-0 3 

3.57925E-03 

3.57950E-03 

3.583 16E-03 

2* 60974E-03 

2.80990E-03 

2.81228E-C3 

2. 10377E-03 

2# X0393E-03 

2. 10554E-03 

1.53472E-C3 

1.53486E-03 

1 • 53600E-C 3 

1.09353E-03 

1 • 09366E-03 

1 • 09446E-C3 


11 111 

N SUB R 
5 


C SUB HX C SUB RHO Y C SUB HY 
• 33333E-01 4. 00000E-Q1 4.25000E-01 

BETA FACTOR 
1.02000E 00 


HL 

SIGMA BAR/E 

HH 

SIGMA BAR/E 

AXISYM 
SIGMA BAR/E 

2.63472E-02 

2.63407E-02 

5.443316-02 

1.562126-02 

1.56155E-C2 

3.842346-02 

9 • 99283E-03 

9.987826-03 

2.721666-02 

6. 840046-03 

6 .835056-03 

1.B9883E-02 

5. 17599E-03 

5. 171576-03 

1 • 34403E-02 

4.180896-03 

4.17611 E— 03 

9.331396-03 

3.58323E-03 

3* 57948E— 0 3 

6.665286-03 

2 • 81231 E-Q3 

2.80990E-03 

4.66569E-03 

2.10557E-03 

2.10392E-03 

3. 29898E-03 

1.536016-03 

1 • 53487E-03 

2. 332856-03 

1.09448E-03 

1. 093656-03 

1 • 63299E-03 
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Figure 11 - Flow-Diagram - Program 4267 
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Figure 11 - Flow- Diagram - Program 4267 
(Continued) 
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Figure 11 - Flow-Diagram- Program 4267 
(Continued) 
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TABLE V - Fortran Listing - Program 4267 


SIBFTC REA00 M94/2,LIST 

COMMON/CSTRAT/DUM ( 50 ) 

C *********************************************************************** 


c standard input for all options 

COMMON/ALL/ 

C 

* TITLE(IO) 

C 

* »RA0RST 
C 

* t NCASES 
C 

*»P0ISR 

C 

* » EUSIG 
C 

* t CALPHL 

c 

*»CALPHU 

c 

♦fCBETAL 

C 

♦fCBETAU 


*»NRRHOY 
*rlETAP 
*»IC12 
*t IC21 
* r 1A12 


*>NR 


*» I DUMP 
♦fORDMAX 


PROBLEM IDENTIFICATION 

INPUT OPTION (RATIO OR STIFF) 

NUMBER OF CASES 

POISSON’S RATIO 

E/SIGMA 

LOWER C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
UPPER C SUB BETA 

THE NUMBER OF (R/RHO SUB Y) RATIOS 
TO BE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEQUE 
NCING OPERATIONS. WILL BE LEFT BL 

ank when only point solutions are 
to be obtained 

ETA SUB P OPTION (BLANK OR 1) 

C SUB 12 OPTION (BLANK OR 1) 

C SUB 21 OPTION (BLANK OR 1) 

A SUB 12 OPTION (BLANK OR 1) 

N SUE R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
MINIMIZATION PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 

SCREENING CUT OFF FOR THE MINIMIZ 


AT I ON PROCESS 


*» SCREEN 


c 


. minimum number of circumferential 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 

HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 

*'CHWMIN 'DUMP 

logical dump 

DUMP=, FALSE. 

DATA RATI /5HRATIO/ 

*»STIF /5HSTIFF/ 

*'BLANK/6H / 

15=5 

16=6 

16=8 

100 CONTINUE 

ASSIGN 1SU0 TO IEOF 
CALL ERCTRL(34,IE0F) 

CALL ERRTR(-34) 

READ(I5'200)TlTLt 
200 FORMAT U0A6) 

CALL ERCTRL(34»0) 

CALL £RRTR(34> 

DUMPS, FALSE. 

RE AD ( 1 5 » 30 Q ) R AORST ' NC ASES » PO I SR r EDS I G ' C ALPHL » C ALPHU r CBETAL 'CBETAUr 
♦NRRHOYr IETAP» IC12» IC2l» IA12,NR» IDUMP 
300 FORMAT (A5»I5'6(F5«0)»7(I5)) 

READ ( 15 » 40 0)ORQMAX» SCREEN »CHWMIN 
400 FORMAT (3E10.5) 

IF ( I DUMP . EQ , 1 ) DUMP= . TRUE • 

IF(RAORST.EQ.RATl >GO TO 600 
IF(RA0RST.EQ.STIF )GO TO 1000 
WRITE(I6»500) RAORST 

500 FORMAT (1HO,45HERROR IN OPTION REQUESTED NOT RATIO OR STIFF »A5) 
RETURN 

600 CONTINUE 
950 CALL RATIO 
GO TO 100 
1000 CONTINUE 

CALL STIFF 
GO TO 100 
1500 RETURN 
END 

SIBMAP TAPt 

ENTRY .UN08. 

•UN08, P2E UNIT06 

UNIT08 FILE * A ( 1 ) 'READY »0UTPUT»8LK=115#bCD 

END 

SIBFTC RATIOD M94/2»LIST 
SUBROUTINE RATIO 
COMMON/ INTER/ 

*F1 

*»F2 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 


* » F3 

*»F4 

*»F5 

*»F6 

*»F7 

*>Fe 

*»F9 

* » FlO 

*»FU 

* » F12 

*»F13 

* » F14 


*»F15 
* » SOLU 1(8) 

*♦ ARRAY <11 '25) 

* * FLHOX( 11 ) 

* » FMMAX 

C*****4****^************************************************************ 

standard input for all options 

COMMON/ ALL/ 


c 

c 

c 


* TITLL(IO) 

* t RAORST 

* t NCASES 
*»POISR 
♦rEDSIG 

* r CALphL 
* » CALPHU 
* » CBLTAL 
♦rCBETAU 


*»nrrhoy 

* t IETAP 
* » IC12 
*f IC21 


problem identification 

INPUT OPTION (RATIO OR STIFF) 
NUMBER OF CASES 

poisson* s Ratio 

E/SIGMA 

LOWeR C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
UPPER C SUB BETA 

THE NUMBER OF (K/RHO SUB Y) RATIOS 
TO aE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEQUE 
NCING OPERATIONS. WILL BE LEFT BL 
ANK WHEN ONLY POINT SOLUTIONS ARE 
TO BE OBTAINED 

ETA SUB P OPTION (BLANK OR 1) 

c sub 12 option (Blank or i> 

C SUB 21 OPTION (BLANK OR i) 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 


C 

*»IA12 


*»NR 


*» I DUMP 
♦rORDMAX 


A SUB 12 OPTION (BLANK OR 1) 

N SUB R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
MINIMIZATION PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 

SCREENING CUT OFF FOR THE MIN1MIZ 
AT I ON PROCESS 


*» SCREEN 


*f CHWrtlN 
LOGICAL DUMP 


MINIMUM NUMBER OF CIRCUMFERENTIAL 
HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 
♦ DUMP 


RATIO INPUT INFORMATION 
COMMON/CSTRAT/ 

CASE NUMBER 


* NCASt 

*» ITYPEO 
♦rTHICKX 

*»flgrad 

♦rTHICKY 

♦♦CRHOX 

*»CHX 

♦»CRHOY 

*»CHY 

*r STRING (2) 


OUTPUT OPTION 1=T ABLEST 2=TABLES 
AND PLOTS , 3=PL0TS» 4=P0INT SOL. 

THICKNESS RATIO (T BAR SUB X/T) 

LENGTH TO RADIUS RATIO (L/R) 

THICKNESS RATIO (T BAR SUB Y/T) 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB RHO X 

GEOMETRIC FACTOR (ONLY FOR POINT 
COLUTION) C SUB H X 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB RHO Y 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB H Y 

STRINGER LOCATION ( INSIDE » OUTSIDE 

♦symmetric 

RING LOCATION ( INSIDE»OUTSIDE»SYM 
METRIC) 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 


**RINGL(2) 

(L/RHO SUB X) LONGITUDINAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 

£ f FLRHOX 

(R/RHO SUB Y) CIRCUMFERENTIAL SLE 
NDERNESS RATIO (POINT SOL. ONLY) 

* 9 frrhoy 

INCREMENTING FACTOR USED IN MINI- 
MIZATION PROCESS 

*»BETAF 

(L/RHO SUB X Y) ARTIFICIAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 

* rFLROXY 

(Cll/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL, ONLY) 

*»Cllk 

(C22/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL, ONLY) 

*»CE2k 

(C12/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL, ONLY) 

* * C 12K 

(C21/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

R/RHO SUB Y RATIOS TO BE IN PLOTS 
OR TABLES 

* » C2lK »RDRH0Y(25) 

DIMENSION CH2(2) 

*»AUT0(2) 

logical point 

DATA ( AUTO ( I ) f 1 = 1 , 2 ) /GHAUTOMA » 3HT IC/ 

♦»FINgID/EHINSIDE/ 

*»0UTSID/6H0UTSID/ 

* * SYMMET/6HSYMMET/ 

P0INT=. FALSE. 

PI=3. 1415926 

PI2=PI*PI 

XCOUNT=0 

IFLAG=0 

15=5 

Ib=6 


itt=8 

FLR0X(2)=‘+2.5 
FLROX ( 3 ) -bO . 0 
FLR0X(4)=66.U 
FLROX (5) =121 • 5 
FLROX (6) = 175, 0 
FLROX(7)=245.0 
FLROX ( B ) =350 . 0 
FLROX (9) =495.0 
FLROX (10)=700.Q 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 


FLROX (11 ) -10o0. 0 
100 CONTINUE 
IENDsll 
FLR0X(1)=30.0 
IF ( NRRHO Y • £0 • 0 ) GO TO 200 
READ (I5» 150) (RDRHOY(I) ,I=1 ,nRRHOY) 

150 FORMAT ( 8E1 0 • 5 ) 

200 CONTINUE 

DO 5000 UU=1'NCASES 

READ ( 15 » 250 ) NCASE » IT YPEO » THICKX » FLGRAD t THICK Y » CRHOX *CH2» CRHOY »CHY 
250 FORMAT ( I5> I5*4ElG,5t A6» A4»2E10.5) 

ICOUNT=ICOUNT+l 
DO 300 I=l»2 

IF(CH2(I) .NE.AUTO(I) )GO TO 350 
300 CONTINUE 

CHX= (1. 0/2. 0)*(THICKX-1.0) /THICKX 
GO TO 450 

350 READ400»NCASE»lTYPEOrTrtICKX»FLGRAD»THlCKY#CRHOX*CHX#CRHOY*CHY 
400 FORMAT ( I5» I5»7E10.5) 

450 CONTINUE 

READ (15 r 50 0 ) STRING » R INGL » BETAF 
500 FORMAT ( A6»A4»A6»A4»Ei0.5) 

if(itypeo.ne.4>go to 1000 

RE AD (1 5 1 550 ) FLRHOX » FRRHOY » FLROXY , Cl 1R t C22R r C 1 2R » C2 1R 
550 FORMAT ( 7E10.5) 

rdrhoy=frrhoy 

900 P0INT=. TRUE. 

1000 CONTINUE 

IF(ITYPEO.eQ,3) GO TO 1050 
CALL PRINT (IFLAG) 

1050 CONTINUE 

DO 2000 I=l»NRRHOY 

F1=PI2/ (RDRHOY ( I ) *RDRHOY ( I ) *FLGRAD*FLGRAD) 

IF(DUMP)WRITE(I8»10000)F1 

F2=0,0 

IF(POINT)F2=(2,0*PI2)/(EDSIG*(THICKX**.5>*(THICKY**.5)*FLROXY^FLRO 

*XY) 

IF(I£TAP.Eq.o)F 2=0.0 
IF(DUMP)WRiTE(I8»l0l00)F2 
F4= ( 4 . 0*PI 2 ) / (FLGR AD*FLGRAD ) 

IF (DUMP) ttRlTE( I8» 10300 )F4 

F 5= ( 1 . 0/ ( 2 . 0*PI2 ) ) *FLGR AD*FLGR AD* ( THICKY** . 5/THICKX** • 5 ) 

IF (DUMP) *RlTE( 18# 10400 )F5 
IF(STRING.EQ.FINSID)FK3X=1.0 
IF (POINT ) GO TO 1100 
IF(STRING.EG.0UTSID)FK3X=-1.Q 
IF ( STRING .£G .SYMMET ) FK3X=0 . 0 
IF(RINGL.EQ.FINSID)FK3Y=1.0 
IF(RINGL.E0.OUTSID)FK3Y=-1.0 
IF ( RINGL .EQ . SYMMET ) FK3Y=0 . 0 
C22R= (FK3Y*CHY ) / ( CRHOY*RDRHOY U ) ) 
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1100 


1200 


1900 

2000 

5000 

10000 

10100 

10200 

10300 

10400 

10500 

10600 


TABLE V - Fortran Listing - Program 4267 
TCont inued) 


C12R=0.0 

C21H=0.0 

CONTINUE 

F7=< , 5*C12r* ( THICKY**, 5/THICKX** • 5 ) ) 

IF(IC12.El,0)F7=0.0 

IF(DUMP)WRlTE<l8,l0600)F7 

F8=( ,b*C2lR*(THlCKX**,5/THICKY**.5) ) 

IF(IC2l.EQ,0)Fd=0.0 
IF (DUMP) WRITE (18 r 10700 )F8 

F9=(-1.0)*FOI5R*( (ThICKY**. 5) / (THICKX**. 5) ) 
IF(IAi2.EL,0)F9=0.0 
IF(DUMP) WRITE (18 » 10600 )F9 
F10=EDSIG* (THICKX**. 5)* (THICKY**. 5) 

IF (DUMP) WRITE (18, 10900) F10 

F12= ( 1 . 0/P 12 ) *FLGRAD*FLGRAD* ( THICK Y/THICKX ) 

IF (DUMP ) WRITE (18, UlOO )F12 
F13=PI2*(C12K/(FLGRAD*FLGRAD) ) 

IF( ICl2.EU.0)Fl3=0.0 

IF ( DUMP ) WRITE (18,1 1200 )F13 

F15=( (PI2*PI2)/(FLGRAD**A.0) )*C12R»C12R 

IF (DUMP) WRITE (18, 11400 )F15 

IF ( , NOT .POINT ) GO TO 1200 

IEND-1 

FLKOXll)=FLRHOX 
uo i9uo u = i , i E nd 
F3=PI2/(FLR0X(J)*FLR0X( J) ) 

IF (DUMP ) WRITE (13,10200 ) F3 

IF ( . NOT . PO I NT ) Cl 1R= ( (FK3X*CHX ) /CRHOX ) *FLGRAD* ( 1 , O/FLROX ( U ) ) 
F6= .b* (CHR+C22R) 

I F ( DUMP ) wR I TE ( 1 8 , 1 0 50 0 ) F6 
F11=PI2 /(FlR0X(J)*FLK0X(U) ) 

IF ( DUMP ) WRITE (18, HOOO)F11 

F14= (PI2*PI2*THICKX ) / (FLGRAD*FLGRAD*THICKY*FLROX l J ) *FLROX ( J) ) 

IF(DUMP) WRITE ( 18, 11300 )F14 

IOPT=0 


CALL MINMUM(IOPT) 

I F ( I T YPEO . EQ . 2 . OR . ITYPEO . EQ . 3 ) ARRAY ( J » I ) =SOLUT ( 3 ) 
IF< ITYPEO. EQ. 3) GO TO 1900 
CALL PRINT2(J»I) 


CONTINUE 

CONTINUE 


IF ( ITYPEO. 
CONTINUE 
RETURN 
FORMAT (1H0 
FORMAT (1H0 
FORMAT (1H0 
FORMAT (1H0 
FORMAT (1H0 
FORMAT ( 1H0 

FORMAT (1H0 


EQ.2.0R. ITYPEO. EO, 3 

, 10X,6HF1 

= » 1PE12 

,1oX»6HF2 

= , 1PE12 

,10X,6HF3 

= t 1PE12 

,10X,6HF4 

= * 1PE12 

,1GX,6HF5 

= » 1PE12 

,10X,6HF6 

r , 1PE12 

,10X,6HF7 

= t 1PE12 


)CALL 


.5) 

• 5) 
.5) 
.5) 

• 5 ) 

• 5) 
.5) 


PLOT 
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TABLE V - Fortran Listing - Program 4267 
(Continued) 

10700 FORMAT tlHO f 10X»6HF8 = ,1PE12.5> 

10800 FORMAT (1HO,1oX»6HF9 = '1PE12.5) 

10900 FORMAT(lH0 v l0X»6HFl0 = »1PE12.5) 

11000 FORMAT ( 1H0 » 1QX » 6HF 1 1 = »1PE12.5) 

11100 FORMAT <1H0,1QX»6HF12 = »1PE12.5> 

11200 FORMAT<lHO'10X'6HFl3 = »1PE12.5) 

11300 FORMAT (1HO,10X,6HFU = * 1PE12.5) 

11400 FORMAT* 1HO,10X»6HF15 = *1PE12.5) 

END 

SIBFTC STIFFD M94/2,LIST 
SUBROUTINE STIFF 

c*:m>******************************************************************** 

C INPUT FOR STIFF OPTION 
COMMON/CSTRAT/ 

C 

* NCASE 
C 

*» All 
C 

»»A22 
C 

**A12 
C 

*»A33 
C 

*»C11 
C 

*rC22 
C 
C 

*»R 
C 


c 

*»D11 

ELASTIC CONSTANT 

c 

*»D22 

ELASTIC CONSTANT 

c 

*»D12 

ELASTIC CONSTANT 

c 

*»D33 

ECCENTRICITY COUPLING CONSTANT 

c 

*»C12 

ECCENTRICITY COUPLING CONSTANT 

c 

*rC21 

OVERALL LENGTH OF CYLINDER 

c 

*»FLENGT 

AN INCREMENTING FACTOR USED IN MI 

c 


NIMIZATION PROCESS 


*»betaf 

• 


COMMON/ INTER/ 
*F1 



CASE NUMBER 

ELASTIC CONSTANT 

ELASTIC CONSTANT 

ELASTIC CONSTANT 

ELASTIC CONSTANT 

ECCENTRICITY COUPLING CONSTANT 

ECCENTRICITY COUPLING CONSTANT 

RADIUS TO MIDDLE SURFACE OF BASIC 
CYLINDRICAL SKIN 

ELASTIC CONSTANT 
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TABLE V — Fortran Listing — Prograa 4267 
(Continued) 


*»F3 

*»F4 

*»F5 

* t Fb 
* » F7 
*»F8 
*»F9 
* » Flu 

* »F11 

* r F12 
* » F13 
*»F14 
*»F15 

* » SOLUT ( 8 ) 
*»ARRAY(llr25) 

* > FLROX (ID 

* »FMMAX 


C ******************** *************************************************** 

standard input for all options 

COMMON/ ALL/ 


* TITLL(IO) 

* t kaorst 

* t NCASLS 
* » POISR 
* t EDSiO 
*»calphl 

* » CALPHU 
*»CbtTAL 
♦ » CdtT AU 


PROBLEM IDENTIFICATION 

INPUT OPTION (RATIO OR STIFF) 

NUMBER OF CASES 

POISSON’S RATIO 

E/SIGMA 

LOWER C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
UPPER C SUB BETA 


THE NUMBER OF (R/RHO SUB Y) RATIOS 
TO BE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEGUE 
NCING OPERATIONS. WILL BE LEFT BL 
ANK WHEN ONLY POINT SOLUTIONS ARE 
TO BE OBTAINED 

*»NRRHOY 

C ETA SUB P OPTION (BLANK 

* t IETAP 

C C SUB 12 OPTION (BLANK 

* * IC12 

C C SUB 21 OPTION (BLANK 


OR 1) 
OR 1) 
OR 1) 


6-30 

GENERAL DYNAMICS CONVAIR DIVISION 



o o o o o o o no non 


GDC-DDG-67-006 

TABUS V - Fortran Listing - Program 4267 
(Continued) 


*» IC2i 
C 

IA12 


*»NR 


♦rIDUMP 
• t ORDMAX 

SCREENING Cl3T OFF FOR THE MINIMI2 
AT I ON PROCESS 

*r SCREEN 

MINIMUM NUMBER OF CIRCUMFERENTIAL 
HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 

* »CHWMIN »DUMP 

LOGICAL DUMP 
15=5 
16=6 
18=8 
IFLAG=0 
100 CONTINUE 

DO 1000 K=1#NCASES 

READ(l5'150)NCASE'Ali,A22’Al2'A33'Cll'C22rR 
150 FORMAT(H0,7E10.5) 

READ(I5»2Q0)DU»D22#Di 2»D33»C12»C21»FLENGT 
200 FORMAT (10X#7ElO#5) 

READ (15 #250) BET AF 
250 FORMAT (10X»E10,5) 

PI=3. 1415926 

FLENG2=FLENGT*FLENGT 

PI2=PI*PI 

R2=R*R 

Fl= ( PI2*D22*A22 ) /FLENG2 

IF(IC12.EQ.0)F1=0.0 

IF (DUMP) WRITE(I8#1000Q)F1 

F2= ( 2 , 0*PI2* ( D12+2 . 0*D33 ) * ( ( A1 1* A22 ) ** . 5 ) ) /FLENG2 

IFUETAP.Eq.Q)F2=0.0 

IF(DUMP) WRITE (I8»lOlOO)F2 

F3= ( PI2*DU*A11 ) /FLENG2 

IF (DUMP) WRITE ( 18 » 10200 ) F3 

F4=(4.Q*PI2)/ (FLENG2/R2) 

IF (DUMP) WRITE (18# 10300 )F4 

F5= (1.0/(2.0*PI2) )*(FLENG2/R2)* ( ( AU/A22)**.5) 
IF(DUMP)WRITE(I8,10400)F5 
F6=(1.0/2.0)*((C11/R)+(C22/R) ) 

IF (DUMP) WRITE ( 18# 10500 )F6 


A SUB 12 OPTION (BLANK OR 1) 

N SUB R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
MINIMIZATION PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 
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F7=(1.0/2.0)*(C12/R)*( ( A11/A22) ** .5) 

IF(IC12.EG.0)F7=0.0 

1F(DUMP) WRITE (18, 10600 )F7 

F6=( 1.0/2. o)*(C2l/R)*( ( A22/A11 ) ** .5) 

IF(IC2l.EO.Q)F8=0.0 

IF ( DUMP) WRITE (18, 10700 )F8 

F9=(2.0*Ai2)/( (All*A22)**.5) 

AF(IAl2.EG,0)F9=0.0 
IF(DUMP) WRITE ( 16, 108U0)F9 
FlO=A33/( (All*A22)**.5) 

IF (DUMP) WRITE (16, 10900) F10 
Fil= (PI2*DU*A11 ) /FLENG2 
IF (DUMP) WRITE ( I 8, 11 000) FI 1 
F12=(1.0/F'12)*( FLENG2/R2 ) * ( A11/A22) 

IF ( DUMP ) WRITE (18, lll00)F12 
F13=(PI2*C12*R)/FLENU2 
IF ( DUMP ) WRITE ( 18, 11200 )F13 
Fl4=(Pl2*Pl2*R2*Dll*A22)/(FLENG2*FLENG2) 

IF ( DUMP ) WR ITE ( I 8 , 1 1300 ) F14 
F16=(PI2*FI2*R2*C12*C12)/(FLENG2*FLENG2) 

IF(DUMP) WRITE (18, ll4U0)F15 
IOPT=l 

CALL MINMUM (I OPT) 

CALL PRINT 3 (K ) 


1000 CONTINUE 
RETURN 


10000 

FORMAT 

(1H0 

,lUX,6HFl 

z 

* 1PE12 • 5 ) 

10100 

FORMAT 

(1H0 

, 10X»6rtF"2 

r 

* 1PE12.5) 

10200 

format 

(1H0 

, 1 OX , 6HF 3 

z 

* 1PE12.5) 

10300 

format 

(IMG 

, 10X , 6HF4 

z 

, 1PE12.5) 

10400 

format 

( 1 HU 

,lUX,6HF5 

z 

* 1PE12 . 5 ) 

10500 

format 

( 1 HO 

,10X,6HF6 

z 

* 1PE12 . 5 ) 

lObOO 

format 

(1HU 

, 1(JX»6HF7 

z 

t 1PE12.5) 

10700 

format 

(1H0 

,10X,6HF8 

z 

» 1PE12.5) 

1060C 

format 

(1H0 

,10X,6HF9 

z 

* 1PE12 . 5 ) 

10900 

F ORMAT 

(1H0 

,1UX,6HF10 

z 

t 1PE12.5) 

11000 

format 

(1H0 

,lOX,6HFll 

z 

t 1PE12 • 5 ) 

111QU 

format 

( 1H0 

,10X,6HF12 

z 

t 1PE12 .5) 

11200 

format 

(1H0 

, 10X, 6HF13 

z 

» 1PE12 • 5 ) 

11300 

format 

(1H0 

, 10 X , 6HF I 4 

z 

* 1PE12.5) 

11400 

format 

end 

(1HQ 

, lUX , 6 HFI 5 

• 

» 1PE12. 5) 

5IBFTI 

C MINIMD N 94/2, LIST 




SUBROUTINE MINMUM (I OPT) 

C***** ********** ************** ************** **************************** 

C STANDARD INPUT FOR ALL OPTIONS 

COMMON/ ALL/ 

C PROBLEM IDENTIFICATION 

* TITLE(IO) 

C INPUT OPTION (RATIO OR STIFF) 

* » Ragrst 
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C 

C 

C 

C 

C 

C 

C 


* fNCASES 

*»P0ISR 

*»EDSIG 

*»CALPHL 

*»CALPHU 

♦fCBETAL 

*»CBETAU 


* » NRRHOY 
*t IETAP 
*f IC12 
* t IC21 
*t IA12 

*rNR 

*» I DUMP 
♦rORDMAX 

*t SCREEN 


*»CHWMIN 
LOGICAL dump 
COMMON/ INTER/ 


*F1 
*»F2 
*»F 3 


NUMBER OF CASES 
POISSON’S RATIO 
E/SIGMA 

LOWER C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
UPPER C SUB BETA 

THE NUMBER OF (R/RHO SUB Y) RATIOS 
TO BE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEQUE 
NCING OPERATIONS. WILL BE LEFT BL 
ANK WHEN ONLY POINT SOLUTIONS ARE 
TO BE OBTAINED 

ETA SUB P OPTION (BLANK OR 1) 

C SUB 12 OPTION (BLANK OR 1) 

C SUB 21 OPTION (BLANK OR 1) 

A SUB 12 OPTION (BLANK OR 1) 

N SUB R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
MINIMIZATION PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 

SCREENING CUT OFF FOR THE MINIMIZ 
ATION PROCESS 

MINIMUM NUMBER OF CIRCUMFERENTIAL 
HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 
.DUMP 
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*»F4 

*»F5 

*>F6 

* rF7 

* t Fd 

*tF9 

* > Flu 
* » Fll 
*rFl2 

*tFXi 

*»F 14 
*»F15 
*»S0LUT(8) 

* f ARRAY (11 >25) 

*»FLkGX(11> 

* t FMMAX 

C********* ****** ******************************************************** 

C RATIO INPUT INFORMATION 

COMMuN/CSTRAT/ 

C CASE NUMBER 

* NCASE 

OUTPUT OPTION 1=TABL£S» 2=TABLES 
AND PLOTS » 3=PL0TS» 4=P0INT SOL. 


*f ITYPEO 

* t THICKX 
*»FLGRmD 
* » T HICKY 

* » CRHOX 

* t ChX 

* t CRHO Y 
* » CHY 

*» STRING ( 2 ) 
*»RINGL(2) 

* » flrhox 


THICKNESS RATIO (T BAR SUB X/T) 

length to radius ratio (l/r> 

THICKNESS RATIO (T BAR SUB Y/T ) 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB RHO X 

GEOMETRIC FACTOR (ONLY FOR POiNT 
COLUTION) C SUB H X 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB RHO Y 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB H Y 

STRINGER LOCATION ( INSIDE»0UTSIDE 
r SYMMETRIC 

RING LOCATION ( iNSlDErOUTSIDEfSYM 
METRIC) 

(L/RHO SUB X) LONGITUDINAL SLENDE 
RNEbS RATIO (POINT SOL. ONLY) 
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(R/RHO SUB Y) CIRCUMFERENTIAL SLE 
NDERNESS RATIO (POINT SOL. ONLY) 

*»FRRH0Y 

INCREMENTING -FACTOR USED IN MINI- 
MIZATION PROCESS 

**BETAF 

(L/RHO SUB X Y) ARTIFICIAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 

•rFLROXY 

(CU/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

*»C11R 

(C22/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

*»C22R 

(C12/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

* »C12R 

(C21/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

R/RHO SUB Y RATIOS TO BE IN PLOTS 
OR TABLES 

*»C21R »RDRH0Y(25) 

EQUIVALENCE 

♦ ( ITYPEO# All ) 

*' (THICKX»A22) 

*t (FLGRAD»A12) 

*t ( THICKY » A33) 

*t (CRHOX.Cll) 

*» (CHX.C22) 

*» (CRHOY »R) 

*» (CHY»D11) 

*» (STRING(I) »022) 

*» (STRING(2) »D12) 

*» (RINGL(i) rD33) 

*» (RINGL(2)»C12) 

*» (FLRH0X»C21) 

* » ( FRRHO Y t FLtNGT ) 

DIMENSION BET ( 9 ) » YMP ( 9 ) » FMMP ( 9 ) 

18=8 
KOUNT=0 
PI=3. 1415926 
FMMIN=1.0 

IF (DUMP) WRITE (I8rl0000)FMMlN 

FMMAX= ( SCREEN/F3 ) ** . b 

IF(FMMAX.LT.2.0)FMMAX=2.0 

IF(QUMP)WRlTE(I8»lOlOO)FMMAX 

6ETMIN=(F1/SCREEN)**.25 

IF (DUMP) WRITE (18. 10200 )BETMIN 

IF(IOPT)30o»300»400 

300 BETAUL=(FMMAX*(THICKX*#.25) )/(FLGRAD*.5*CHWMIN*(THICKY**.25) )*PI 
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60 TO 450 

400 BETAUL=(FMmAX*PI*R*(A 22**.25) )/<FLENGT*.5*CHwMlN*(All** .25) ) 

450 FMAXSY = ( 1 • 0/ (F14+F15) ) ♦*•25 
IF (DUMP ) WHITE (18, 10300 )BETAUL 
I F ( FMAXSY . LT . 1 . 0 ) FMAXS Y=1 . 0 
IF (DUMP) WRITE (I8» 10400 )FMAXSY 

YAXSYW=F11*FMAXSY*FMAXSY+ (F12/ (FMAXSY*FMAXSY ) ) * ( ( 1 . Q-F13*FMAXSY*FM 
♦AXSY ) **2 • 0 ) 

IF ( DUMP )WRITE( 18 rl 0500 )YAXSYM 
IF ( I OPT )50o#500»600 
500 S0LUT(8)=YAXSYM 
60 TO 700 

600 S0LUT(8)=YAXSYM/AU 
700 CONTINUE 

bETA=BETAUL 
1100 bETA2=BETA*B£TA 

tlETA4=BETA2*BETA2 

IF (DUMP) WRITE (18, 10600) BETA 

bETAS=BETA 

CALL 5EVELE(BETAS»YS»FMS»I0PT) 

K0UNT=K0UKT+1 
IF(K0UNT.Lt.l)60 TO 1650 
IF(YS.6E.YS1) GO TO 1700 
1650 YS1=YS 

IF ( DUMP ) WHITE (18, 11000)YS1 
FMSI=FMS 

IF (DUMP) WRITE (18, 111U0)FMSI 
bETASl=BETAS 

IF(DUMP) WRITE ( 18, 1 1200) BETAS1 
1700 IF (BETAS .LE.BETMIN)G0 TO 1800 
bETA=BET A/BET AF 
60 TO 1100 
1600 CONTINUE 

bETAFl=BETAF 
DO 2100 I=1»NK 
bETAF2=SQRT (BETAF1 ) 

IF (DUMP) WRITE ( 18, 113U0)BETAF2 

DO 19U0 J=l»9 

FJ=-5+U 

bET( J)=BETAS1/BETAF2**FJ 
IF(BET( J) .LE.BETAUL)60 TO 1850 
BET(J)=BETAUL 

1650 CALL SEVELE(BET(J) »YMP( J) ,FMMP(U) ,IOPT) 

1900 CONTINUE 
K=1 

DO 2000 J=l,9 
IF(YMP(J) .LT.YMP(K) )K=J 
2000 CONTINUE 

YS1=YMP(K) 

IF(DUMP) WRITE ( 18,1 1400) YS1 
BETAS1=BET(K) 
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IF ( DUMP > WR ITE ( 1 8 , 1 1500 ) BET ASl 
FMSI=FMMP(k) 

IF(DUMP)WRlTE(I8»Il6U0)FMSI 
BETAF1=BETAF2 
2100 CONTINUE 

IF( IOPT)22oO» 2200 *2500 
2200 S0LUT(3)=YS1 
SOLUT U»=FMSI 

SOLUT ( 2 > = < FMSI*PI* ( THICKX** , 25 ) ) / ( BET AS1*FLGRAD* < THICKY** .25 ) ) 
IF ( YSl.GT • YAXSYM)G0 to 2300 
60 TO 2400 

2300 S0LUT(3)=YAXSYM 
SOLUT(l)=FMAXSY 
SOLUT (2) =0 , 0 
2400 GO TO 3000 
2500 CONTINUE 

IF(YS1.GT.YAXSYM)G0 TO 2800 

S0LUT(3)=YS1/A11 

SOLUT(l)=FMSI 

SOLUT (2 > = < FMS I *PI *R* ( A22** . 25 ) ) / ( BET AS1 *FLENGT* ( A1 1** .25 ) ) 

60 TO 3000 

2800 SOLUT (3>=YaXSYM/AU 
SOLUT (1)=FMAXSY 
SOLUT (2) =0.0 

3000 IF(YS1.GT.YAXSYM)G0 TO 3500 
FML=CALPHL*FMSI 

fmh=calphu*fmsi 
bETL=CBETAL*BETASl 
8ETH=CBET AU*BET AS1 
SOLUT ( 4 ) =FCN ( BETL » FML ) 

SOLUT ( 5 ) =FCN < BETH » FML ) 

SOLUT ( 6 ) =FCN ( BETL » FMH ) 

SOLUT ( 7 ) =FCN ( BETH » FMH ) 

IF ( IOPT)35oO» 3500 » 3100 
3100 00 3400 1=4 » 7 

SOLUT ( I ) =SOLUT ( I ) /All 
3400 CONTINUE 
3500 RETURN 

10000 FORMAT(1HO,10X'12HM SUB MIN = * 1PE12.5) 

10100 FORMAT (1H0,10X»12HM SUB MAX = ,1PE12.5) 

10200 FORMAT(lHO,loX»15HBETA SUB MIN = »1PE12.5) 

10300 FORMAT(1HO,10X»14HBETA SUB UL = ,1PE12.5) 

10400 FORMAT (1HO,10X»15HM SUB AXISYM = #1P£12.5) 

10500 FORMAT C1H0,10X»15HY SUB AXISYM = »1PE12.5> 

10600 FORMAT(1HO,10X» 7HBETA = »1PE12.5> 

11000 FORMAT(lHO,lOX,2oHlTERATlVE Y SAVED = rlPE12.5) 

11100 FORMAT ( 1H0 » 1 OX » 2 OH ITERATIVE M SAVED s »1PE12.5) 

11200 FORMAT ( 1H0 f 1 OX *23H ITERATIVE BETA SAVED = >1PE12.5) 

11300 FORMAT(1HO,10X»20HBETA FACTOR SUB 2 = »1PE12.5) 

11400 FORMAT (1HQ,10X»31HREFINEMENT ITERATION SAVED Y = »1PE12.5> 
11500 FORMAT (1HO,1oX# 34HREFINEMENT ITERATION SAVED BETA = rlPE12.5) 


6-37 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG- 6 7-006 


TABLE V - Fortran Listing - Program 4267 
(Continued) 


llfaOO FORMAT <lHO,10Xr31HREFINEMENT ITERATION SAVED M = ,1PE12.5) 

END 

SIBFTC SEVELD *94/2, LIST 

SUBROUTINE SEVELE(B£TAS»YS»FMStIOFT) 

COMMON/ INTER/ 

*Fl 
*» F2 
*»F3 

* r F4 

* ,F5 
*,F6 

* »F7 
* , F6 
*»F9 
*»Flu 
*»F11 

**FU 

* > F13 
*»F14 
* » Fib 

* t SOlu r ( 8 ) 

* , Array ( H * zb ) 

* r FLkOa (11) 

*»FMMAx 

£********************************************************♦************** 

C RATIO INPUT INFORMATION 

COMMOH/CS I RAT/ 


c 

c 

♦ NCAbt 

CASE NUMBER 

OUTPUT OPTION 1=TABLES» 

2=TABLES 

c 


AND PLOTS f 

3=PL0TS» 4=P0INT SOL. 

c 

*» ITYPEO 

THICKNESS 

RATIO (T BAR 

SUB X/T) 

c 

* , THICRX 

length TO 

RADIUS RATIO 

(L/R) 

c 

*f F lORAD 

THICKNESS 

RATIO (T BAR 

SUB Y/T) 

c 

* f ThIOKY 

GEOMETRIC 

FACTOR (ONLY 

FOR POINT 

c 

C 

*r CRhOX 

SOLUTION) 

GEOMETRIC 

C SUB RHO X 
FACTOR (ONLY 

FOR POINT 

c 

c 

*, CHX 

COLUTION) 

GEOMETRIC 

C SUB H X 
FACTOR (ONLY 

FOR POINT 

c 

C 

* t CRHOY 

SOLUTION) 

GEOMETRIC 

C SUB RHO Y 
FACTOR (ONLY 

FOR POINT 

c 

*>CHY 

SOLUTION) 

C SUB H Y 
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♦♦STRING (2) 

*»RINGL(2) 

♦♦FLRHOX 


STRINGER LOCATION ( INSIDE ♦OUTSIDE 
♦SYMMETRIC 

RING LOCATION ( INS IDE ♦ OUTSIDE ♦SYM 
METRIC) 

(L/RHO SUB X) LONGITUDINAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 

(R/RHO SUB Y) CIRCUMFERENTIAL SLE 
NDERNESS RATIO (POINT SOL. ONLY) 


♦♦FRRHOY 

INCREMENTING FACTOR USED IN MINI- 
MIZATION PROCESS 

♦♦BETAF 

• (L/RHO SUB X Y) ARTIFICIAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 

♦♦FLRQXY 

(Cll/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 


♦♦CiiR 


(C22/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 


♦♦C22R 

(C12/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL, ONLY) 

♦♦C12K 

(C21/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

R/RHO SUB Y RATIOS TO BE IN PLOTS 
OR TABLES 

♦♦C21R rRDRH0Y(25) 

^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦** 
C STANDARD INPUT FOR ALL OPTIONS 
COMMON/ ALL/ 

C PROBLEM IDENTIFICATION 


C 

C 

C 

C 

C 

C 

C 


♦ TITlE(IO) 

♦♦RAORST 

♦♦NCASES 

♦♦POISR 

♦♦EDSI6 

♦♦CALPHL 

♦♦CALPHU 

♦♦CBETAL 


INPUT OPTION (RATIO OR STIFF) 
NUMBER OF CASES 
POISSON»S RATIO 
E/SIGMA 

LOWER C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
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C 

*»CBETAU 


*»nrrhqy 

* » IETAP 
*> IC12 
* * IC21 
*» IA12 


>MNR 


* » I DUMP 
♦rORDMAX 


*» SCREEN 


♦rCHWMlN 
LOGICAL DUMP 
EQUIVALENCE 
* ( ITYPEO* A1 1 ) 

*» (THICKXr A22) 

*» (FLGRAD»A12) 

*» (THICKY # A33) 

*t (CRHQX»Cll) 

*r (CHX»C22> 

*» (CRHOY »R> 

*» (CHYrDll) 

*> ( STRING ( 1 ) #022) 
*» (STR1NG(2) »012> 
*» (RINbL(l) »D33) 
*» (RINGLU) f Cl2) 
*» (FLRH0X»C21) 

*t (FRRHOY»FlENGT) 
PI=3. 1415926 


UPPER C SUB BETA 

THE NUMBER OF (R/RHO SUB Y) RATIOS 
TO BE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEGUE 
NCING OPERATIONS. WILL BE LEFT BL 
ANK WHEN ONLY POINT SOLUTIONS ARE 
TO BE OBTAINED 

ETA SUB P OPTION (BLANK OR 1) 

C SUB 12 OPTION (BLANK OR 1) 

C SUB 21 OPTION (BLANK OR 1) 

A SUB 12 OPTION (BLANK OR 1) 

N SUB R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
MINIMIZATION PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 

SCREENING CUT OFF FOR THE MINIMIZ 
AT I ON PROCESS 

MINIMUM NUMBER OF CIRCUMFERENTIAL 
HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 
rDUMP 
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PI2=PI*PI 

BETA=B£TAS 
BETA2=BETA*BETA 
BETA4=BETA2*BETA2 
IF ( IOPT ) 80 o , 800*900 

800 FMSl=(BETA*FLeRAD*CHWMIN*CTHlCKY**t25) )/(PI*2.0*(THICKX**.25) ) 

GO TO 1000 

900 FMS1=<BETA*FlENGT*CHwMIN*(AU**.25> )/ (PI*R*2.0*(A22«‘*.25) ) 

1000 CONTINUE 

IF(FMSl.LT.1.0)FMSl=1.0 

IF(DUMP)WRITE(I8,10700)FMS1 

0MIC=.5*< (F4/Fl)**.5)*(F6*(1.0/(BETA2 ) >+F7+F6*(1.0/(BETA4 ) 

* ) ) 

IF (DUMP) WRITE (18, 10800 )0MIC 

FM=( ( (F5**2)*F4*4.0*BETA4*BETA4)/(4.*F1*< (1.0+(F2/F1)*BETA2+(F3/F1 
* ) *BET A4 ) * ( 1 . o+F9*BET A2+F1 0*BET A2+BET A4 ) +4 . 0*OMIC*OMIC*BETA4*BETA4 ) 

IF (DUMP )WRlTE(l8r 10900 )FM 
IF(FMS1.LT,FM.AND,FM.LT.FMMAX>G0 TO 1200 
Y1=FCN(BETA»FMS1) 

IF (DUMP) WRITE ( 18, 11000)Y1 
Y3=FCN( BETA* FMMAX) 

IF (DUMP) WRITE ( I8,Hl00)Y3 
Y 2= ABS ( Y1 ) ABS ( Y3 ) 

GO TO 1300 

1200 Y1=FCN(BETA»FMS1) 

IF(DUMP)WRiTE(I8,11000)Y1 

Y2=FCN(BETA*FM) 

IF (DUMP) WRITE( 18, 11200 )Y2 
Y3=FCN ( BET A * FMMAX ) 

IF (DUMP) WRlT£(I8,lllO0)Y3 
1300 CONTINUE 

IF(Y1.LT.Y2.AND.Y1.LT.Y3)G0 TO 1500 

IF(Y2.LT.Y3)G0 TO 1400 

FMS= FMMAX 

bETAS=BETA 

YS=Y3 

GO TO 1600 
1400 FMS=FM 

betas=betm 

YS=Y2 

GO TO 1600 
1500 FMS=FMS1 

bETAS=BETA 

YS=Y1 

1600 CONTINUE 

IF (DUMP) WRITE (18, 11300) BETAS 
IF (DUMP) WRITE (18,1 1400 )FMS 
IF ( DUMP ) WRITE (18,1 1500 )YS 
RETURN 
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10700 FORMAT ( 1H0 » l0X» 18H (M SUB S) SUB 1 = tlPE12,5) 

10800 FORMAT ( 1H0 , 10X» lOHUPSILON = »1PE12.5) 

10900 FORMAT (lHQ»10Xr 4HM = »1PE12.5) 

11000 F0RMATUHG,1QX,37HY CORRESPONDING TO (M SUB S) SUB 1 = #lPE12.b> 
11100 FORMAT ( 1H0 » loXf 31HY CORRESPONDING TO M SUB MAX = »1PE12.5) 

11200 FORMAT(lHO,lOXf23HY CORRESPONDING TO M = flPEl2.t>) 

11300 FORMAT(1HO,10X,13HSAVED BETA = #1PE12.5) 

11400 FORMAT <1H0,10X,1QHSAVED M s »1PE12.5) 

11&00 FORMAT(1HO,10X»10HSAVED Y r »1PE12,5) 

END 

5 IBFTC Y M94/2»LIST 

FUNCTION FCN(BETA»FM) 

COMMON/INTER/ 

*F1 

*>F2 

*»F3 

*>Fh 

* » F5 

*rF6 

* » F7 

*»F8 

*»F9 

*»F10 

*»F11 

*fFl2 

*rFl3 

*»F14 

*»Flb 

* t SOLUT ( 8 ) 

* r ARRAY (H'25) 

*»FLROX(ll) 

* rFMMAX 
18=8 

bETA2=BETA*BETA 
bETA4=BET A2*8ETA2 

FCN=Fi*( (FM*FM)/BETA4)+F2*( (FM*FM ) /BETA2) +F3* (FM*FM> +F4*FM*FM*8ETA 
*4* ( F5* ( 1 • / ( FM*FM ) ) -F6* ( 1 . /BETA2 ) -F7-F8* < 1 . Q/BETA4 ) ) **2 / ( 1 . Q+F9*B 

*ETA2+F10*BETA2+BETA4 ) 

RETURN 

END 

SItiFTC PRINTD M94/2»LIST 

SUBROUTINE PRINT (IFLAG) 

COMMON/INTER/ 

*F1 

*»F2 

*»F3 

*»F4 

* t Fb 
*»F6 


6-42 

GENERAL DYNAMICS CONVAIR DIVISION 



o o no on on no on on on non on 


GDC-DDG-6 7-006 


TABLE V - Fortran Listing - Program 4267 
(Continued} 


*rF7 

* *F8 

* *F9 

* *F10 
*»FU 
*r Fl2 

* »F13 
*»F14 
*»F15 
*#S0LUT(8> 
*»ARRaY(X1»25) 


*»flrox(U> 

*»FMMAX 

C*********************************************************»************* 
C KATIO INPUT INFORMATION 

common/cstrat/ 

c CASE NUMBER 

* NCASE 

OUTPUT OPTION 1=TABLES» 2=TA6LES 
AND PLOTS* 3=PL0TS» 4=P0INT SOL. 


*r ITYPEO 
♦rTHICKX 

* »FLGRaD 
*»THICKY 

* r CKHGX 
*rCHX 
♦rCRHOY 

* rCHY 

*»STRInG(2) 

*»RIN6L(2> 


THICKNESS RATIO <T BAR SUB X/T ) 

LEN6TH TO RADIUS RATIO (L/R) 

THICKNESS RATIO (T BAR SUB Y/T) 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB RHO X 

GEOMETRIC FACTOR (ONLY FOR POINT 
CQLUTION) C SUB H X 

geometric factor (only for point 

SOLUTION) c SUB RHO Y 

GEOMETRIC FACTOR (ONLY FOR POINT 
SOLUTION) C SUB H Y 

STRINGER LOCATION ( INSIDEfOUTSIDE 
* SYMMETRIC 

RING LOCATION (INSIDE»0UTSIDE»SYM 
METRIC) 

(L/RHO SUB X) LONGITUDINAL SLENDE 
RNESS RATIO (POINT SOL. ONLY) 


* »FLRHOX 

(R/RHO SUB Y> CIRCUMFERENTIAL SLE 
NDERNESS RATIO (POINT SOL. ONLY) 

**FRRHOY 
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*»BETAF 

*»flrqxy 

*»cuk 

*»C22R 
* » C12K 


*»C21R 


incrementing factor used in mini- 
mization PROCESS 


(L/rtHO SUB XT) ARTIFICIAL SLLNDE 
RNESS RATIO (POINT SOL. ONLY) 

(Cll/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

(C22/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 


(CI2/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 


(C21/R) ECCENTRICITY COUPLING RAT 
10 (POINT SOL. ONLY) 

R/RHO SUB Y RATIOS TO BE IN PLOTS 
OR TABLES 
»RDRHOY ( 25) 

********************************** 


C 

C 


STANDARD INPUT FOR ALL OPTIONS 
COMMON/ ALL/ 

* TITLE(IO) 

♦rRAORST 

* » NCASES 

* » POISR 

*»EDSIG 

♦rCALPHL 

* » CALPHU 

*»CBETAL 

*»CbETAU 


* »nrrhoy 
ietap 


problem identification 

input OPTION (RATIO OR STIFF) 

NUMBER OF CASES 

POISSON'S RATIO 

E/S1GMA 

LOWER C SUB ALPHA 
UPPER C SUB ALPHA 
LOWER C SUB BETA 
UPPER C SUB BETA 

THE NUMBER OF (R/RHO SUB Y) RATIOS 
TO BE INCLUDED IN PLOTS OR TABLES 
WHICH RESULT FROM AUTOMATIC SEGUE 
NCING OPERATIONS. WILL BE LEFT BL 
ANK WHEN ONLY POINT SOLUTIONS ARE 
TO BE OBTAINED 

ETA SUB P OPTION (BLANK OR 1) 

C SUB 12 OPTION (BLANK OR 1) 
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C 

c 


** IC12 
**IC2l 
* * I A12 

**NR 

*» I DUMP 
*»ORDMAX 


C SUB 21 OPTION (BLANK OR 1) 

A SUB 12 OPTION (BLANK OR 1) 

N SUB R *** THE NUMBER OF REFINE- 
MENT CYCLES TO BE EMPLOYED IN THE 
minimization PROCEDURE 

DUMP OPTION *** 1 WHEN SUPPLEMENT 
ARY DIAGNOSTIC OUTPUT DESIRED 

MAXIMUM ORDINATE FOR PLOTS 

SCREENING CUT OFF FOR THE MINIMIZ 
AT ION PROCESS 


** SCREEN 

minimum number of circumferential 

HALF WAVES USED TO ESTABLISH A 
CUT OFF VALUE IN THE MINIMIZATION 
PROCEDURE 

•rCHWMlN >DUMP 

LOGICAL DUMP 
EQUIVALENCE 
*(ITYPEO»All) 

*» (THICKX»A22) 

** (FLGRAD* Aj.2) 

*# (THICKY»A33) 

** (CRHQXfCll) 

*» (CHX»C22> 

*» (CRHOY#R) 

* * (CHY fDll) 

(STRXNG(l) »D22) 

** <STRING(2> *D12> 

*» (RINGL(l) r 033) 

*i (RINGL(2) »Cl2 ) 

*» (FLRH0X»C21) 

*» (FKRHOY»FLENGT) 

I5=S 

16=6 


IF(ITYPE0.EQ.3)RETURN 
IFLAG=IFLAG+1 
IF(IFLAG.Ne.1)g0 to 200 
WRITE(I6»5000> 

WRITE ( 16 »SiOO) TITLE 
WRITE ( 16 >5200) 

WRITE ( 16*5300) 

WRITE ( 16*5400) R AORST * NC ASES * PO I SR * EDS I G * C ALPHL * C ALPHU > CBET AL rCBETA 
*U»NRRHOY* l£TAP» IC12» 1C21* IA12 


WRITE ( 16*5500) 


6-45 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 


TABLE V - Fortran Listing- Program 4267 
(Continued) 


WRITL(I6»b600) 

1F(ITYPE0.£G.4)G0 TO 100 

WRITE. ( 16*5700) ORDMAx * SCREEN , CHWM IN » NR 

WRITE, ( 16 » bQOO ) 

WRITE(I6*59U0) (RDRHOY(I) *I=1*NRRH0Y) 

60 TO 200 

100 CONTINUE 

WRITE ( 16 * 7500) SCREEN *CHWMIN* NR 
200 CONTINUE 

IF ( IFLAG.Ne.1) WRITE ( 16*7750 ) 

WRITE ( 16* 600u ) 

WRITE ( 16* 6100 ) 

IF ( IT yPE0»EG»4) GO To 400 
1F(CH2.EQ.BLANk)G0 TO 300 

WRITE (16* 6200 )NCASEr IT YPEO,THICKX*FLGRAD*THI cKY*CRHOX*CHX»CRHOY*CH 
*Y 

sO TO 50u 

DATA BLANK/6H / 

300 CONTINUE 

WRITE (16* 7600 )NCASE,lTYPEO,THlCKX»FLGRAD»THlCKY»CRHOX»CRHOY»CHY 
GO TO 500 
400 CONTINUE 

WRITE (16 *7700) NCASE*1TYPE0,THICKX*FLGRAD*THICKY 
wRi it (16*7300) 

wR1Te(I6*740u)FLKHOx»FRRHOY»FLROXY*C 11R*C22R»C12R»C21R 
WRITE ( I6*t,3Uu) 

wR ITE (16 *6400) STRING »RINGL»BETAF 
RETURN 

600 wRIl E ( 16 * 630u ) 

wRlTL(I6»6400)STRlNG*RINGL»tiETAF 

RETURN 

5000 FORMAT (1H1,31X,62HG£NERAL INSTABILITY OF ECCENTRICALLY STIFFENED C 
♦YLINUERS 4267) 

5100 FORMAT ( 1H0 t 32X» 1oA6) 

5200 FORMAT ( 1HG , 5HINPUT * 4X » 6HNUMBER * 3X * 8HP0I SSONS „ 13X » 5HL0WER » 8X » 5HUPPE 
*R * 7X f 5HL0v.£R , 7X » BHUpPtR# 13X* 9HETA SUB P *4X r 3HC12* 5X * 3HC21 » 5X * 3hA12 
* ) 

5500 FORMAT (1X»6H0PTI0N*2X,BH0F CASESt 3X*5HRATI0*5X»3HE/G»4X*llHC SUb A 
*LPHAt2X,llHC SUB ALPHA* 2X* 10HC SUB BETA»2X,l0HC SUB BETAr2X*6HNRRH 
*OY,5X*6HOPTION»4X»6HOPTlONr2X*6HOPTIONf2X»6HOPTION) 

5400 F0KMAT(lX*A5,5X*14*5XfF5.3f2X»F7.3»5X*F6.3»7x,F6.3*7x*F6.3*6X»F6.3 
**6X» I2»7 a» I1,9X»I1»7A» I1*7X»11) 

5500 FCRMmT(1Hu,22X,12HMaX ORDINA TE* 10X»9hSCREENING» 15X » HHWIIN NO CIRC) 
5600 FORMAT (25 x,9HF0R PLOTS* 12X, 7HCUT-OFF* 16X* 10HHALF-WAVES* 11X * 7HN SUB 
* R) 

5700 FORMAT ( 1HQ » 22X» IPElO • 3* 11X» 1PE12.5* 15X 1 0PF5»2* 16X * 12) 

5600 FORMAT (XH0,42Xr40HR/RHO Y RATIOS FOR AUTOMATIC SEQUENCING) 

5900 FORMAT (lHo» (6(1PE12.5,4X) »/lX) ) 

6000 FORMAT (lH0,2X»4HCASE*9X*bH0UTPUT) 

6100 FORMAT (2X»6HNUMB£R*8X,6H0PTI0N,7X#9HT BAR X/T »8X»3HL/R*8Xf 9HT BAR 
*Y/T»5 a*UhC SUb RHO a f 4X * 8HC SUB HX*5X»11HC SUB RHO Y»4X»8HC SUB H 
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♦ Y ) 

6200 F0RMAT<1H0,2X»I4,11X»I1»8X,7(1PE12.5,2X)) 

6300 FORMAT (1H0,21X,17HStRINGER LOCATION# 17X,13NRING LOCATION# 16X # 11HBE 
♦TA FACTOR) 

6400 F0RMAT(1H0,25X»A6,A3»25X»A6»A3»17X»1PE12.5) 

6500 FORMAT (1HO,5«X#7HR/RHO Y) 

6600 FORMAT ( 1H0 * 55X, 1PE12* 5) 

6700 FORMAT aH0,l6X»8HCRlTlCAL#6X#8HCRITICAL#6X»8HCRITICAL»8X»2HLL#l2X» 
♦2HLH » 12X » 2 hHL » l2x # 2HHH » 10X # 6HAXI ; SYM ) 

6600 FORMAT (3X#7HL/RHO X,5X, 12HL0WER CASE M#2X» 12HL0WER CASE N,2X,6(UH 

♦SIGMA BAR/E»3X)> 

6900 FORMAT ( 1H0 ,9 ( 1PE12»5»2X) ) 

7000 FORMAT (1H0,9X»6HCRITICAL»6X»8HCRITICAL»6X#8HCRITICAL»8X#2HLL#12X#2 
♦HLH # 12X # 2HHL , 12X # 2HHH # 1 OX # 6HAX I S Y« ) 

7100 FORMAT < 8X »12HL0WER CASE M»2X, 12HL0WER CASE N,2X»6( UHSIGMA BAR/E #3 

♦ X) ) 

7200 FORMAT ( 1H0 * 7X#8( 1PE12 t 5, 2X) ) 

7300 FORMAT ( 1H0 »28X,7HL/RH0 X#7X#7HR/RH0 Y#7X»8HL/RH0 XT#6X , 5HC11/R,9X, 
♦5HC22/R * 9X , 5HC12/R » 9X » 5HC21/R ) 

7400 FORMAT UH0 f 27X#7(lPEl2.5#2X)) 

7500 FORMAT ( 1H0 , 43X # 1PE12 . 5 # 15X # F5 . 2 # 16X » 12 ) 

7600 FORMATClHO f 2X»l4,llX»Il»8X#4(lPEl2*5»2x)»14X#2(lPE12.5#2X)) 

7700 FORMAT (lH0#2Xr 14# 11X» II #8X#3J1PE12.5#2X)) 

C 

^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦* *♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
£♦*♦*♦*♦*♦*♦*♦♦♦*♦*♦*♦*♦♦♦♦♦*♦♦♦*♦♦♦*♦*♦*♦♦♦*♦♦♦*♦*♦♦♦*♦♦♦*♦♦♦♦*♦♦♦♦♦♦♦♦ 

C 

ENTRY PR1NT2 ( JJ, I I ) 

IFUTYPE0.EQ.4) GO TO 700 
IF(JJ.NE.l) GO TO 600 
IF ( I I ,NE . 1 ) WRITE ( 16,7750 ) 

7750 FORMAT (1H1) 

WRITE ( I6#650u) 

WRITE ( 16,6600 IRDRHOYC II) 

WRITE ( 16,6700) 

WR I TE( 16,6800) 

600 CONTINUE 

IF(S0LUT(2).NE.0.0) GO TO 650 

WRITE ( I6»7775)FLR0X (jj) , (SOLUT ( I ) , 1=1 , 3) ,S0LUT(8) 

7775 FORMAT (1H0,4(1PE12»5,2X) , 56X, 1PE12«5) 

RETURN 

650 WRITE(I6»8900)FLROX(Jj) , (SOLUT ( I ) , 1=1, 8) 

RETURN 

700 WRITE (16,7000) 

WRITE(I6,71Q0) 

IF(S0LUT(2) .NE.0.0) Go TO 750 
WRITE(I6,725u) (SOLUT(I) ,1=1,3) ,S0LUT(8) 

7250 F0RMAT(1H0,7X,3(1PE12.5,2X) ,56X, 1PE12.5) 

return 

750 WRITE(I6,720C) (SOLUT( I ) , 1=1,8) 

RETURN 
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C 

c ******************************************* **************************** 
c ************************************* ********************************** 

c 

entry PRINJ3 (KK ) 

ITYPEO=0 

15=5 

16=6 

16=6 

1F(KK.NE.1)G0 TO I0o0 
WRITE ( I6» 5000 ) 

WRITE ( 16 » 5100 ) TITLE 
WRITE ( I6» 7800 ) 

WRITE ( 16 » 7900 ) 

WRITE (I6f 6000 )RA0RST»NCASES»CALPHL»CALPHU»CBETAL»CBETAU»IETAP»IC12 

* t IC21 > IA12 
WRITE ( 16 * &10Q ) 

WRITE ( 16 » 6200 ) 

WRITE ( 16 » 6300) SCREEN rCHWMIN»6£TAF» NR 
1000 WRITE ( 16 » 6400 > 

WRITE ( 16 » 6500 ) 

WRITE ( 16 » 6600 )NCASE,ITYPE0, All »A22»Al2,A33,ClltC22»R 
WRllLll6r8700) 

WRITE(I6»660U)D11»D22,D12»D33»C12»C21»L 
wRITe( I6»690u) 

WRITE ( 16 » 9000 ) 

WRITE(I6»9l00) (SOLUT(I) »I=lf6) 

RETURN 

7600 FORMAT ( 1H0 , SHINPUT » 7X , 6HNUMBER 1 6X » 5HL0WER * 8X , 5HUPPER 1 8X » 5HL0WER » 7 X 
*»5HUPPER»7x»9HETA SUE Pt 4X , 3HC12»9X » 3HC21 t 9X » 3HA12 ) 

7900 FORMAT (1X»6H0PTIQN»5X»8H0F CASES » 2X r 11HC SUB ALPHA »2X» 11HC SUB ALP 
*HA,2X»10HC SUB BETA,2X»1QHC SUB BETA»7X»6H0PTI0N»4X»6H0PTI0N»6X»6H 
*UPT10ivi» 6X > 6H0PTI0N ) 

6000 FORMAT(lHG,A5»8Xrl4,7x»F6.3»7X>F6.3»7X»F6.3»6X,F6.3»lOX*Il»9X»Il»l 

*lX»IlfllX»ll) 

6100 FORMAT(1HO,18X,9HSCREENING,11X» hhMIN NO C IRC » 19X * 4HBET A ) 

8200 FORMAT ( IX » 19X »7HCUT-0FF * 12X# 10HHALF— WAVES 1 19X »6HF ACTOR » 11X» 7HN SUB 

* R ) 

8300 FORMAT UH0,16X,lPEl2.5r 12X»F5.2» 18X» 1PE12.5# 11X» 12) 

8400 FORMAT ( 1H0 , IX » 4HCASE » 3X 1 6H0UTPUT ) 

8500 FORMAT < IX 1 6HNUMBER t 2X , 6H0PTI0N » 6X » 3HA11 1 11X » 3HA22# 11X » 3HA12# 11X f 3H 
*A33»llX»3Hcll»llX»3HC22t 12X»1HR) 

8600 FORMAT (1H0,1X»I4»5X, I 1,5X#7(1PE12.5»2X) > 

8700 FORMAT(lHO,2oX»3HDllrllXr3HD22»llX»3HDl2»llXr3HD33»llX»3HCl2#llX»3 
*HC2l » 12X » lHL ) 

8800 FORMAT ( 1H0,16X» 7 (1PE12.5»2X) ) 

8900 FORMAT (lH0,2X»8HCRlTICAL»7X»8HCRITICALr5X#8HCRITICAL»10X»2HLL#13X» 
*2HLH»l3X»2HHL»l3X»2HHH»10Xr6HAXlSYM) 

9000 FORMAT (IX* 12HLOWER CASE M» 3X» 12HL0WER CASE N»3Xr6(7HN SUB X»8X)) 
9100 F0RMAT(1Hu,1PE12.5»3x,4(1PE12.5#2X) » IX # 3( 1PE12.5» 3X ) ) 

end 
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SIBFTC PLOTD *94/2»LIST 
SUBROUTINE PLOT 
COMMON/ INTER/ 

*F1 
*»F2 
*rF3 , 

*»F4 

*»F5 

*»F6 

* »F7 
*»F8 
*rF9 
*»F10 
*>F11 
*»F12 
*»F13 
*»F14 
*»F15 

*»SOLUT(B) 

* t ARRAY ( 1 1 » 25 ) 

**FLROX<U) 

♦»fmmax 

C************************************************* 1 *** ,,, ***** 1, ‘ ,,t * 1,t ****** 1,tl,t * 

COMMON/ ALL/ 

* TITLEdO) 

*»RAORST 

♦rNCASES 

*»POISR 

♦fEosie 

*»CALPHL 

*»CALPHU 

♦»CBETAL 

♦rCBETAU 

*»nrrhoy 
*t IETAP 
*»IC12 
* » IC21 
*»1A12 
*»NR 
*» I DUMP 
*»ORDMAX 
*» SCREEN 

*»chwmin rDUMP 

logical dump 

c*********************************************************************** 

common/cstrat/ 

* NCASE 
♦rlTYPEO 
* » THICKX 
♦rFLGRAD 
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♦rTHICKY 

♦*CRHOX 

♦ »CHX 

♦ *CRHOY 
* * CHY 

♦*STRING(2) 

♦*RINGL(2) 

♦»flrhox 

♦»FRRHOY 

♦*betaf 

♦ * flroxy 

♦*C11R 

**C22R 

♦ f C 1.2.R 

♦*C21R >RDRH0Y(25) 

t-XTERNAL N XV 
EXTERNAL NYV 
MX=1 
MY=U 

CALL SMXYV(MXrMY) 

L=4 

CALL SETMIV(24»0,130»24) 

XL=10.0 

XR=1000.0 

YB=0.0 

YT=ORDMAX 

UX=1,0 

IF(ORDMAX.LT..U0fa.AND.ORDMAX.GT..001)DY=.0001 

IF(ORUMAX.GE..006)DY=.0002 

IF (ORDMAX.lE. .001) DY=. 00002 

N=1 

IF(DY.EQ. .0001)M=10 

IF(DY.EQ. .0002. OR. DY.EQ. .00002 )M=5 

I=-N 

J— —Ni 

iMX=6 

NY=b 

CALL GRIDlV(L»XL.XR» YB»YT*DX*DY*N*M»I* J*NX,NY) 

CALL APRNTV(0*-14»-21,21HSIGMA BAR / ELASTIC E»5»790) 

CALL PRINTV(-14»14HT BAR X / T = *133*1010) 

CALL LABLV(THlCKX*245, 1010.-4* 1*1) 

CALL PRINTV(-6,6HL/R = .551*1010) 

CALL LABLV(FLGRAD» 599,1010.-4*1*1) 

CALL PRINTV(-14»14HT BAR Y / T = .823*1010) 

CALL LABLVITHICKY *935 .1010,-4*1*1) 

CALL PRINTV (-9.9HL / RHO X, 554*120) 

CALL RlTE2v(248*77*lU23»90.1»38*-l»3dHCRITICAL AXIAL COMPRESSION F 
♦OR etNERAL.NLAST) 

CALL RITE2V (239*45* 1023*90*1* 39 »-l* 39HINST ABILITY IN CYLINDERS WIT 
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*H ECCENTRIC tNLAST) 

CALL RITE2V < 392 # 13 » 1023 » 90 # 1 r 22 » -1 > 22H0RTH0TR0PIC STIFFENING » NlAST 
*1 

00 1000 I=1#NRRH0Y 

CALL APL0TV(11»FLR0X»ARRAY(1»I)»1#1»1»55»IERR) 

1000 CONTINUE 
RETURN 
ENO 

SIBMAP ERRTR1 
* 

* 

* USAGE- FORTRAN Iv 

♦ • 

♦ • 

♦ CALL ERRT* (ICODE) 


♦ # 

* « 

♦ • 

♦ WHERE 

♦ 

♦ 

* 

* REFERENCE 

♦ 


IABSUCODE) = FORTRAN EXECUTION ERROR MONITOR CODE 
IF CODE .lT. 0 * THEN NO ERROR TRACE WANTED. 

IF CODE .GT. 0 » THEN ERROR TRACE WANTED. 

* COMPUTER SYSTEMS BULLETIN » NO. 90 » 9 MAY 66 


(FXEM) 


ERRTR 

SAVE 

2.4 


AXT 

U.2 


CLA* 

3.4 


ST2 

flag 


TPL 

again 


STL 

flag 


SSP 


AGAIN 

SUB 

=35 


TXI 

♦+1.2.-1 


TZE 

*+2 


TPL 

again 


STA 

als 


CLA 

=1 

ALS 

als 

** 


nzt 

FLAG 


TRA 

off 


ORS 

ErtRw-1 »2 


RETURN 

ERRTR 

OFF 

COM 



ANS 

ERR#— 1 #2 


RETURN 

ERRTR 

•ERRW. 

CONTRL 

ErRW»ERRWNo 

ERRW 

PZE 



PZE 



PZE 


ERKWND 

NULL 
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* 

* 

FLAG 


TABLE V - Fortran Listing - Progran 4267 
(Continued) 

IF FLAG ,EQ. 0 * then turn BIT OFF. 

IF FLAG ,NE. 0 t THEN TURN BIT ON, 

BSS 1 
END 
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APPENDIX A 

BUCKLING CURVES FOR SANDWICH CYLINDERS 


In order to provide a basis for the comparison of skin-stringer- 
ring constructions against sandwich configurations, the buckling curves 
of Figure 12 are presented here. These curves are based on the classical 
solution of Stein and Mayers C 17] for axially compressed, simply supported 
sandwich cylinders having isotropic faces of equal thickness and an 
isotropic core. Since this solution employs small— deflection shell theory, 
the values obtained from these curves must be reduced by suitable knock- 
down factors to account for the influences of initial imperfections (see 
Volume V [l2}) e This practice is necessary if one is to obtain safe design 
values for the critical loading. The following definitions apply to the 
notation used in Figure 12: 
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Critical running axial compression, lbs/in. 
Overall length of cylinder, in. 

Poisson's ratio of faces. 

Mean radius of cylinder, in. 

Distance between middle surfaces of faces, in 
Young's modulus of faces, psi. 

Thickness of single face, in. 

Shear modulus of core, psi. 

Core thickness, in. 


A- 2 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 



Figure 12(a) - Classical Buckling Coefficients For Simply 
Supported Sandwich Cylinders Under Axial 
Compression (Isotropic Core And Faces) 


A- 3 

GENERAL DYNAMICS CONVAIR DIVISION 







GDC-DDG-67-006 



z 

a 

Figure 12(b) - Classical Buckling Coefficients For Simply 
Supported Sandwich Cylinders Under Axial 
Compression (Isotropic Core And Faces ) 
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APPENDIX B 

BUCKLING CURVES FOR ISOTROPIC CYLINDERS 

In order to provide a basis for the comparison of skin-stringer-ring 
constructions against isotropic cylinders, the buckling curves of Figures 13 
through 18 are presented here. Each of the following loading conditions 
are covered: 

(a) Unpressurized under pure axial load. 

(b) Unpressurized under pure bending moment. 

These curves were developed by Convair through the statistical analysis of 
test data [l8]. Separate families are given for each of the following 
statistical criteria: 

(a) Best fit 

(b) 90% probability; 95% confidence 

(c) 99% probability; 95% confidence 

Each family consists of separate curves for (L/R) = 0.25, 1.0, and 4.0. 
Criterion (a) was established by means of the conventional least squares 
technique and gives the mean expected (50% probability) level for buckling. 
Half of any large array of typical test data would be expected to fall below 
this level. Criterion (b) may be considered to represent that there is 
95% confidence that at least 90% of any large array of test points would 
fall on or above the related design curves. Criterion (c) may be similarly 
expressed for 99% probability and 95% confidence. Criteria (b) and (c) 
statistically correspond to the MIL-HDBK-5 "B" and "A M values, respectively. 
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Figure 13 - q /E vs. R/t for Unpressurized Monocoque 
Circular Cylinders (Clamped Ends) Under 
Pure Axial Load; BEST FIT 
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R/t 

Figure 14 - g /E vs. R/t for Unuressurized Monocoque 
Circular Cylinders (Clamped Ends) Under 
Pure Axial Load; PROBABILITY = 90 % 
CONFIDENCE s 95% 
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Figure 15 - c^/E vs. R/t for Unpressurized Monocoque 
Circular Cylinders (Clamped Ends) Under 
Pure Axial Load: PROBABILITY ■ 99% 
CONFIDENCE = 95% 
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Figure 16 - /E vs. R/t for Unpressurized Monocoque 
Circular Cylinders (Clamped Enri«) Under 
Pure Bending Moment: BEST FIT 
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Figure 17 - o^/E vs. R/t for Unpressurized Monocoque 
Circular Cylinders (Clamped Ends) Under 
Pure Bending Moment: PROBABILITY = 90% 
CONFIDENCE = 95% 
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Figure 18 - c /E vs. R/t for Unpressurized Monocoque 

C I* 

Circular Cylinders (Clamped Ends) Under 
Pure Bending Moment; PROBABILITY - 99% 
CONFIDENCE 5 95% 
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